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ABSTRACT 


The  purpose  of  the  Lopez  Island  Ocean  Bottom  Seismometer 
Intercomparison  Experiment  was  to  determine  the  effects  of  coupling  and 
bottom  currents  on  ocean  bottom  seismometers.  Twelve  operational  OBS's, 
three  specially  designed  three-component  systems  and  a  hydrophone  were 
compared  with  each  other.  Unlike  seismometers  placed  on  hard  rock  at 
land  stations,  ocean  bottom  seismometers  can  be  affected  by  soft  sediments 
(which  act  as  lossy  mechanical  springs)  and  by  buoyancy.  Coupling  through 
soft  sediments  can  modify  the  response  to  ground  motion  much  as  a  low  pass 
filter  does,  and  high  buoyancy  tends  to  counteract  this  effect.  These 
effects  are  observed  in  the  Lopez  data,  which  consist  of  signals  from 
mechanical  transient  tests,  cap  shots,  airgun  pulses,  and  general  back¬ 
ground  noise.  The  modification  of  response  is  pronounced  for  some 
instruments  and  barely  noticeable  in  others.  Instruments  that  stand  high 
in  the  water  relative  to  their  base  width  tend  to  be  susceptible  to  rock¬ 
ing  motion  that  shows  up  as  a  mechanical  cross  coupling  between  horizontal 
and  vertical  motion.  Correlation  of  Lopez  results  with  coupling  theory 
suggests  that  it  is  possible  to  design  ocean  bottom  seismometers  that  will 
couple  well  to  any  sediment.  Current  levels  at  the  Lopez  site  (<5  cm/sec) 
were  too  small  to  produce  noticeable  effect  on  anv  of  the  instruments; 
however,  the  same  design  criteria  that  will  minimize  coupling  problems  will 
also  lessen  problems  caused  by  ocean  currents. 
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I.  INTRODUCTION 


Between  June  13  and  30,  1978  a  field  intercomparison  of  several 
ocean  bottom  seismograph  (OBS)  systems  was  conducted  in  Shoal  Bay,  Lopez 
Island  (Puget  Sound),  Washington.  Participants  in  the  experiment  are 
listed  in  Appendix  A.  The  two  main  reasons  for  the  field  test  were: 

(1)  to  compare  the  responses  of  the  instruments,  which  have  quite  varied 
mechanical  configurations,  when  coupled  to  the  earth  through  soft  sedi¬ 
ments  comparable  to  those  of  the  ocean  floor;  and  (2)  to  determine  the 
susceptibility  of  the  instruments  to  noise  induced  by  near-bottom  ocean 
currents.  The  results  of  this  experiment  should  provide  information 
important  for  the  design  of  OBS  systems  that  will  minimize  noise  and 
optimize  coupling  characteristics. 

Signals  recorded  on  ocean  bottom  seismograph  geophones  are  often 
monochromatic  and  complicated.  Whether  these  features  are  natural  in 
origin  or  caused  by  the  presence  of  the  OBS  is  a  major  question  to  be 
answered  by  this  test.  Some  data  indicated  that  the  presence  of  the  OBS 
on  the  soft  sediments  could  severely  distort  the  motion  of  the  ground; 
however,  the  extend  of  this  problem  for  OBS's  of  different  configurations 
was  unknown.  The  complications  present  on  many  geophone  records  are  not 
observed  on  ocean  bottom  hydrophones,  thus  the  question:  do  geophone 
data  add  any  useful  information  that  cannot  be  obtained  from  hydrophones? 

Some  data  have  indicated  that  near-bottom  ocean  currents  can  increase 
noise  levels  on  OBS's  enough  to  make  them  unusable  during  high-current 
periods;  therefore,  we  wanted  to  test  as  many  systems  as  possible  for 
susceptibility  to  ocean  current  noise. 

In  addition  to  12  operational  OBS's  supplied  by  10  different  research 
organizations,  three  sets  of  three-component  (T0  =  1  sec)  seismometers 
and  a  hydrophone  were  included  to  provide  standards  for  comparison:  a 
"spike  standard"  was  pushed  firmly  into  the  bottom;  a  "plate  standard" 
with  a  smooth;  hemispherical  superstructure  on  a  large,  flat,  circular 
plate  was  placed  on  the  sediments;  and  a  "neutral  density  standard"  with 
a  roughly  spherical  shape  was  floated  within  the  uppermost  sediment 
(Fig.  1-1).  The  instruments  were  placed  within  a  few  meters  of  each 
other,  and  four  current  meters  were  installed  around  the  array  to  moni¬ 
tor  water  circulation  in  the  bay. 

Data  were  obtained  from  each  sensor  by  hard  wiring  analog  signals 
to  land  where  they  were  recorded  by  a  24-channel  digital  system.  About 
250  seismic  records  were  obtained;  these  include  signals  from  mechanical 
transient  tests,  seismic  signals  from  211  airgun  shots  and  8  blasting 
caps,  and  12  samples  of  background  noise  (Appendix  B).  In  addition  to 
the  digital  recording,  many  of  the  instruments  recorded  internally, 
supplying  valuable  supplementary  data. 

Mechanical  transient  tests  were  conducted  to  provide  some 
quantitative  estimate  of  the  coupling  function  between  the  ocean  bottom 
and  the  OBS  package,  including  possible  cross  coupling  between  horizontal 


and  vertical  motions.  The  procedure  is  analogous  to  the  classical  weight 
lift  test  for  earthquake  seismographs,  although  because  of  buoyancy  the 
response  of  the  seismic  system  to  the  mechanical  transient  test  is  not 
necessarily  the  same  as  that  to  a  seismic  input. 

In  this  report  we  summarize  the  results  of  the  Lopez  Island 
intercomparison  test.  Representative  data  are  shown  with  the  text,  and 
selected  records  from  each  recorded  sensor  are  shown  in  Appendix  C. 
Additional  justifications  for  the  results  are  presented  in  the  referenced 
papers,  included  as  appendices  to  this  report. 


II.  EXPERIMENT  DESCRIPTION 


Achievement  of  the  test  goals  required:  a  logistically  favorable 
site  that  was  reasonably  seismically  quiet  and  had  sediments  with  physical 
properties  comparable  to  those  of  deep-sea  sediments;  a  system  for 
digitally  recording  multiple  channels  simultaneously;  a  method  for 
measuring  total  instrument  response,  including  the  effects  of  the  bottom; 
and  the  measurement  of  physical  oceanographic  data  such  as  tidal  height 
and  current  velocities.  The  implementation  of  each  of  these  requirements 
is  described  below. 


Site  Choice  and  Description 

The  requirements  that  a  site  have  sediments  with  physical  properties 
similar  to  those  of  deep-sea  sediments,  reasonably  low  noise  levels,  and 
be  logistically  accessible  with  respect  to  scuba  divers,  power,  and 
transportation,  removed  from  consideration  most  sites  in  or  adjacent  to 
open  ocean,  because  of  surf  noise,  and  sites  that  were  too  remote.  A 
site  that  met  all  the  requirements  was  situated  at  the  northern  end  of 
Lopez  Island,  in  the  San  Juan  Islands  in  Washington  State.  Its  general 
location  is  shown  in  Figure  II-l  and  the  specific  location.  Shoal  Bay, 
Figure  II-2. 

Shoal  Bay  lies  between  Humphrey  Head  and  Upright  Head,  which 
consist  mostly  of  consolidated  siltstone-sandstone  conglomerates  of 
marine  origin.  The  water  depth  in  Shoal  Bay  is  a  fairly  uniform 
6  to  7  meters  at  low  water  and  is  underlain  by  a  generally  constant  3  m 
of  acoustically  transparent  soft  silty  i..ud;  this  latter  fact  was 
determined  from  reflection  profiles  taken  with  a  3-KHz  system  (Fig.  II-3) . 
It  is  not  known  whether  this  mud  is  underlain  by  the  hard  conglomerates 
or  by  less  competent  glacial  material  that  would  in  turn  be  underlain  by 
the  conglomerates. 

Several  short  gravity  cores  were  taken  in  the  area  of  the 
instrument  array  and  the  physical  properties  of  the  top  7  cm  were  measured 
(Table  1 1- 1) .  These  sediments  were  almost  exclusively  silt  and  were 


V  c 


n  map  of  Lopez  Island  area.  Dash-dot  line  is 
border . 


3.  Reflection  profiles  from  Shoal  Bay.  The  profiles  show  a  relatively  uniform  layer  of  uncon 
lated  sediment  over  a  harder  unknown  material  . 
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TABLE  II-l.  Physical  Properties  of  Shoal  Bay  Sediments 


DENSITY,  g/cc 
POROSITY,  % 


sample:  1 _ 2 _ 3 _ a  ve  rage 

1.60  1.56  1.54  1.57 

71.66  68.03  64.40  68.03 


GRAIN  SIZE 

%  sand  4.3 

7.  silt  9  3.6 
7.  clay  2.0 

MINERALOGY 


Quart  z 

55 

P lagioc  lase 

35 

Mica-Kaolinite- 

10 

heavy  minerals 

found  to  be  slightly  greater  in  density  and  lower  in  porosity  than  typical 
deep-sea  muds. 


Recording  System 

The  general  locations  of  the  OBS  array,  current  meters,  and  the 
house  used  as  recording  base  and  workshop  are  shown  in  Figure  II-2.  The 
relative  locations  of  the  instruments  shown  in  Figure  II-4  and  Table  II-2 
were  initially  determined  by  sea  surface  triangulation  with  respect  to 
the  current  meters  and  later  refined  with  seismic  wave  triangulation 
(Tuthi.ll,  1980).  Instruments  with  only  one  horizontal  geophone  were 
oriented  to  make  the  horizontal  sensitive  to  North-South  motion.  To 
transfer  the  data  from  the  OBS's  to  the  recording  facility  in  the  house, 
a  24-channel  cable  with  underwater  pluggable  connectors  was  deployed. 

Each  OBS  was  required  to  have  matching  connectors  and  cables  that  could 
be  plugged  into  this  cable.  The  shore-based  recording  facility  consisted 
of  analog  recorders  and  a  24-channcl  digital  recording  unit. 


Measurement  of  Response 


In  addition  to  the  recording  of  numerous  noise  samples,  two  types 
of  controlled  source  experiments  were  performed. 

The  first  was  from  number  8  caps  detonated  on  the  bottom  at  distances 
less  than  100  m  from  the  array  and  211  40-cu-in. airgun  shots  fired  at  a 
variety  of  distances  and  azimuths  from  the  array.  These  sources  were 
used  primarily  to  evaluate  the  response  of  the  different  OBS  designs  to 
the  same  source.  Accurate  source-receiver  distances  were  not  determined; 
in  hindsight  this  proved  to  be  an  unfortunate  omission  since  some 
interesting  surface  wave  modes  were  observed  and  dispersion  measurements 
were  degraded  by  the  lack  of  accurate  distance  measurements.  (Informa¬ 
tion  on  source  times  and  locations  is  on  file  with  groups  at  UW  and  HIG.) 

The  second  consisted  of  mechanical  transient  tests  applied  directly 
to  each  instrument.  The  purpose  of  these  transient  tests  was  to  provide 
a  quantitative  estimate  of  the  transfer  function  between  the  ocean  bottom 
and  the  OBS  package,  including  possible  cross  coupling  between  horizontal 
and  vertical  motion.  The  transients  were  obtained  by  electromagnetically 
releasing  a  float  from  each  system,  thereby  generating  a  step  of  force. 
Tests  were  performed  in  the  vertical  and  two  orthogonal  horizontal 
directions  as  shown  schematically  in  Figure  II-5.  Signal  size  was 
adjusted  for  each  instrument  by  modifying  the  number  of  small  wooden 
floats  used  or  the  amount  of  scuba  air  trapped  within  an  inverted  tomato 
can  used  as  a  float.  The  latter  procedure  was  used  for  the  larger 
instruments . 
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.  II-4.  Location  of  OBS  array.  The  array  was  located  by  least-squares  methods  from 
seismic  arrival  time  data.  The  letters  next  to  each  instrument  refer  to  the  insti¬ 
tute  that  built  it  (Appendix  A). 


.  I  I  —  3 .  Mechanical  transient  tests.  Impulses  are  applied  to  each  instrument 
to  determine  coupling  parameters.  A  float  is  released  by  breaking  the  cir¬ 
cuit  in  an  electromagnet  holding  the  float  to  the  instrument.  From  Sutton 
et  a  1  . ,  1 9 MO  (Appendix  F) . 


Physical  Oceanography  Pat  a 

Physical  oceanography  information  gathered  during  the  intercalibration 
study  consisted  of  current  records  from  four  Anderaa  current  meters 
surrounding  the  main  instrument  grouping.  Currents  were  measured  at  a 
distance  of  2  m  above  the  sea  floor  and  recorded  at  15— min  intervals. 

Ihe  current  meter  moorings  were  approximately  40  m  irom  each  other.  Each 
mooring  was  individually  anchored  wiLh  a  single  railroad  wheel  and 
nuoyeci  nv  a  single  ORE  float  oi  positive  buoyancy  (300  lb).  The  currents 
were  extremely  weak;  peak  currents  were  b  cm/sec.  Water  temperatures 
ranged  between  11  and  12  C  and  salinites  ranged  between  35.9  and  36.5  ppt. 

The  S1)S  (sediment  dynamics  system)  was  deployed  near  the  center  of 
the  array.  This  device  records  tide  elevations,  temperature,  and  current 
speeds  and  directions,  and  photographs  the  bottom.  Pressure,  tempera¬ 
ture,  and  currents  were  recorded  at  15-min  intervals.  Photographs  were 
taken  at  1-hr  intervals.  A  limited  number  of  records  are  available  of 
1024  pressure  measurements  taken  it  1-sec  intervals.  Data  are  available 
tor  thi  period  June  20  to  28,  1978.  Because  of  the  weak  currents,  current 
direction  measuri ments  made  with  the  80S  are  not  reliable.  The  speeds 
agree  with  those  measured  with  the  Anderaa  current  meters. 

Tidal  records  tor  the  experiment  duration  are  shown  in  Figure  11-6. 


ill.  I MSTKOIENTATION 


Twelve  instrument  systems  were  tested  at  Lopez  Island  and  their  data 
were  compared  with  data  from  four  "standard"  instruments.  The  characteris¬ 
tics  of  tile  12  instruments  and  tiie  method  used  to  record  data  are 
discussed  in  this  section. 

Each  of  the  12  instruments  records  data  differently  when  in  use  on 
the  ocean  bottom;  however,  since  the  purpose  of  this  test  was  not  to 
judge  recording  systems,  signals  were  taken  from  each  system  as  analog 
voltages  before  the'  normal  recording  systems.  Eacli  institute  supplied 
an  external  plug  on  its  DBS  for  each  sensor,  and  voltages  were  condi¬ 
tioned  to  he  between  - 5  volts.  Wires  from  each  sensor  were  connected 
to  a  junction  box  on  the  bottom  tiiat  was  wired  to  shore.  An  attenuator 
panel  on  land  allowed  for  gain  adjustment  before  digitizing.  Data  were 
digitized  with  a  lb-hit  DFS-IY  multi-channel  system  supplied  by  SI0  and 
were  also  recorded  on  monitor  paper  records.  For  distribution,  the  data 
were  translated  into  a  standard  format  on  9-track  digital  tape.  It  was 
possible  to  record  only  22  sensors  out  of  36  at  any  one  time;  thus  the 
recording  channels  were  rearranged  occasionally  so  that  representative 
signals  could  he  obtained  from  each  component. 

The  standard  instruments  (plate,  spike,  and  neutral  density)  were 
devised  bv  HIC  and  LDGO  to  provide  the  best  fidelity  of  ground  motion 
possible  on  the  basis  of  three  different  assumptions.  The  plate  (PL) 


provides  a  la rp,i. ■  bearing,  surface  and  verv  low  bearing  pressure,  thus 
causing  coupling  resonances  to  be  above  the  seismic  range.  Its  low 
profile  and  streamlined  shape  should  also  minimize  current-induced  noise. 

The  spike  (SI1)  instruments  weredriven  into  the  sediments  to  provide  cou¬ 
pling,  with  material  with  higher  shear  strength  than  that  of  surface 
sediments.  This  should  also  move  coupling  resonances  to  high  frequencies 
and  eliminate  current- induced  noise.  The  neutral  density  standard  (ND) 
was  made  with  a  density  approx  imntelv  that  of  the  sediment  in  which  it 
w.is  buried.  I'll  i  s  instrument  should  respond  without  coupling  effects 
because  ol  the  lack  of  den.xitv  contrast,  and  should  he  insensitive  to 
current- induced  noise  because  ol  its  burial.  Drawings  of  these  instru¬ 
ments  are  shown  in  Figures  111-1  to  111-).  Each  of  the  standards 
contained  three  orthogonal  HSIO-I  geophones  connected  to  a  preamplifier 
(Fig.  ill-Al.  Flectrical  calibration  transient's  were  applied  to  all 
standard  components  to  insure-  proper  (matched)  operation. 

In  addition  to  the  geophone  standards,  a  hvdrophono  standard  was 
provided  In  l.iK.t*  (lie .  Ml-  '.  it  was  supi  1  ied  with  an  internal 
preampl it  i--r . 

The  1 svstcmx  tested  were  as  dill cr«  nt  as  any  set  ol  systems  can  be. 
Drawings  of  each  system  are  shewn  in  alphabetical  order  in  Figures  III  6  to 
17,  and  important  physical  parameters  each  svste::;  are  given  in  Table  III-l 
Blanks  in  the  ta  le  reflect  information  not  available  at  the  time  of  printing 
If  additional  :  n  I  e  r.ia  t  i  on  ivno  >  r.inc  a  pa<*t  tcular  instrument  is  required, 

the  reader  i  s  d  ; t  ,  .  t  .1  :  o  :  r.  M  ■  ;  a  \p,n-tu!  i\  A  g  !••• ;  ng  the  names  and 

addresses  ol  euc  Mist  1 1  HUa  r-.  r-suv.  i:v.  olved  !  n  the  experiment  and 
abbrev  ial  ion  o-  t  :  i .  ■  •  i  tut  'on  and  it.  package. 

Figure  IM-!S  a-ot.-s  the  .relit  ..le  r«  sjum-u  of  each  instrument  used 
in  tin  sv.stem  i  :•>.  •  Ig.li  n-  -eir ,  el-i  Monies,  ar.d  amplifiers  up  to  the 

point  wile r«  the  si., uni  wa-  taken  -'it  lor  the  test.  Although  we  attempted 

to  reference  all  i  ns  1 1  ucieii  t  s  to  the  '-irk  le't  1,  tin  reference  level  may 
he  in  error  in  senu  eases.  Note  also  Mint  tin.-  IIK.  i  espouse  curves  are 

lor  ga  i  n  7  o!  thei*'  automat  ic  i;a  i  o  cent  nd  svsl  iti;  the  absolute  level  i  or 

anv  particular  sistnal  depends  on  the  -.tain  step  in  which  the  sensor  ampli¬ 
fier  was  oper.il  in.-,.  Also,  the  l-n  i  vers  i  t  \  of  Washington  instruments  employ 
square  root  roopres.-,  i  on .  The  signal  -passed  to  the  recorder  is  tin*  square 

root  ol  the  values  shown  in  I  -  air--  Ill-Id.  In  some  cases  where  all 

sensors  have  the  same  response  in  a  p  .it icalar  instrument,  only  one 
response  curve  i  •  -bcon. 

Table  MI-.’  gives  the  w  loci  tv  transti!  function  of  each  instrument 
together  with  the  values  of  associatid  constants. 


pling  with  t ho  sediment .  Four  «.■  1  ongatcd  holes  along  the  outer  edge  of  the  plate  aid 
hindling  .md  emplacement  of  1:k'  instrument.  A  level  bubble  is  emplaced  at  the  top  of 
dome.  Numerous  holes  were  drilled  in  the  plate  and  the  dome  to  allow  any  trapped 
or  or  air  to  escape.  When  in  position  .  t  the  bottom,  the  mudlevel  is  just  below  the 


Fig.  III-2 . 
Spike  Standards 


Two  horizontal  and  one  vertical 
HS-10  1Hz  geophones  were  configured 
as  shown  in  the  sketches.  The  hori 
zontals  were  oriented  N-S  and  E-W. 
All  were  pushed  into  the  sediment 
a  few  CM  below  the  mudline.  The 
spike  frames  were  made  of  1/2" 
aluminum  plate. 


Weight  of  vertical  in  air  =  13.6Kg 

Displacement  =  4.750CC 

Weight  in  sediment  of  pi. 57  =  6.1Kg 


Weight  of  horizontal  in  air  =  11.4K] 

Displacement  *  4,950CC 

Weight  in  sediment  of  pi. 57  =  3.6Kg 
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Standard  hydrophone  deployment. 


Geophysics  Standard  OBS 


Wt .  in  air  -  250  kg. 

Wt .  in  water  -  100  kg 
Bearing  area  -  220  cm*- 
Geophone  HS10-2  located  15  cm  from  end  of 


instrument  cylinder 


Fig.  II'i-9.  Lamont-Doherty  Geological  Observatory  OBS. 


Wt.  In  Air  60  Kg 
Wt.  in  Water  14  Kg 
Ceophone  4.5  Hz 


Fig.  TIT-15.  University  of  Texas  at  Galveston  OBS 


mudline' 


'  W.isliin^tun  '  l-lower  Pot"  OBS . 


RESPONSE,  dB  re  I  volt 


<xl 


i 
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Response  to  Vertical  Input  on  Vertical  Geophones 

The  vertical  motion  of  an  QBS  package  (Fig.  1V-1)  resulting  from  a 
vertical  seismic  input  can  be  approximated  by  the  transfer  function  given 
Equation  1  (Appendix  E) . 


I  =  (l-C)fs2+2h  w  S/(l-C)+w2/(l-C) 

L  C  C  c 

s  — 


2  2 
S  +2h  (i)  S+(u 
c  c  c 


(1) 


where 


1  =  Zj+X^=  displacement  of  instrument  package 
Xj’8  equilibrium  value  of  X 

2^=  displacement  of  water-sediment  interface 
a)  =  angular  frequency  of  OBS-bottom  coupling 


hc=  =  damping  coefficient  of  OBS-bottom  coupling 
C  =  coupling  constant 


Thus,  given  values  for  C,  hc,  and  .c,  the  response  can  be  determined. 
Response  curves  for  three  values  of  coupling  are  shown  in  Figure  IV-2. 

The  coupling  constant  is  given  by 

c  =  (Mi-VV/M*I  (2) 

where 

=  instrument  mass 

M  =  mass  of  displaced  water  at  equilibrium 
w 

M  =  mass  of  displaced  sediment  at  equilibrium 

M  =  "  +M  +M  ,  mass  of  OBS  plus  water  and  sediment  that  move  with  it  . 

I  1  w  s 

The  value  for  C  can  vary  from  zero  at  neutral  bouyancy  where 
=  1  and  seismic  coupling  is  perfect  for  all  frequencies,  to  some 
Slue  less  than  one.  A  maximum  value  for  C,  independent  of  shape,  can 
be  calculated  from 

C  =  1- (M  +M  )/Mt  (3) 

max  w  s  '  1 


=  1/12 


Fig.  IV-2.  Amplitude  response  of  an  OBS  package  (motion  of  OBS/motion  of  bottom)  for 
different  values  of  C.  As  C  approaches  zero,  the  coupling  approaches  optimum.  From 
Sutton  et  al.,  1980  (Appendix  E) . 


(For  most  present  OBS  designs  the  package  is  primarily  in  the 
water  and  Ms  0.)  Values  of  Cmax  calculated  for  the  instruments  compared 
at  Lopez  Island  are  given  in  Table  IV-1.  The  more  flat  the  horizontal  sur¬ 
face  area,  the  greater  C  is  overestimated  by  equation  (3). 

Experimental  values  for  hc  (or  equivalently,  Qc)  and  a>c  can  be 
obtained  from  mechanical  transient  tests  such  as  those  described  in 
section  Il-C.  In  addition,  Mj,  needed  to  estimate  C,  can  be  obtained 
from  the  ratio  of  applied  force  to  output  amplitude.  Equations  (13)  and 
(14)  in  Appendix  K  give  the  transforms  of  the  output  voltage  from  a  criti¬ 
cally  damped  geophone  that  would  be  expected  from  such  tests.  Theoretical 
curves  are  matched  to  Fourier  transforms  of  transient  outputs  from  three 
DBS's  and  the  spike  standard  in  Figure  IV-3.  (A  complete  set  of  transient 
test  outputs  and  spectra  for  each  OBS  and  the  standard  instruments  is  in 
Appendix  C.)  Although  the  observations  are  degraded  by  imperfect  mechani¬ 
cal  testing  procedures,  we  obtained  fairly  reliable  estimates  of  hc  and 
i. c  from  these  comparisons.  Estimates  of  the  coupling  frequency  obtained 
from  the  vertical  transient  tests  are  listed  in  Table  IV-1. 

Theoretical  values  for  Qc  (Qc=l/2hc)  and  mc  can  be  obtained  from 
Figure  1  in  Appendix  E  and  the  data  in  Table  1 1 1  —  2 .  Calculated  mass 
factors,  Fm,  which  provide  an  estimate  of  Qc,  and  calculated  coupling 
frequencies  are  listed  in  Table  IV-1.  These  values  are  based  on  a  single 
estimate  of  sediment  rigidity  and  density,  do  not  include  any  entrained 
sediment  or  water  (M*,  M*=0).  neglect  the  effect  of  buoyancy  on  u'c.  and  do 
not  take  into  account  large  differences  in  the  configuration  of  the  areas 
in  contact  with  the  bottom  (only  the  horizontal  component  of  the  contact 
area  is  considered). 

Observed  and  calculated  coupling  frequencies  are  compared  in 
Figure  IV-4.  Most  of  the  operational  DBS's  show  fairly  good  agreement  for 
an  assumed  sediment  shear  velocity  Vs=10m/s  and  density  ps=l/57g/cm3 . 

LUGO,  UW ,  and  UWF,  however,  fall  well  off  the  line  with  a  higher  observed 
than  calculated  frequency.  From  Table  1 1 1  —  I  we  see  that  these  packages 
apply  a  bearing  pressure  to  the  bottom  sediment  more  than  three  times 
greater  than  that  of  the  next  greatest.  Therefore,  it  is  likely  that 
they  are  coupled  to  a  stiffer  bottom;  e.g.,  as  shown  in  Figure  IV-4, an 
increase  in  Vs  to  33  m/s  would  bring  theory  and  observation  together  for 
LUGO;  smaller  increases  are  required  for  UW  and  UWF. 

All  three  of  the  standard  instruments  show  considerably  lower 
observed  than  calculated  coupling  frequencies.  All  three  instruments  are 
tightlv  coupled  to  the  bottom  sediment  and  it  is  likely  that  there  is 
considerable  error  in  neglecting  M*.  As  indicated  in  Figure  IV-4,  calcu¬ 
lated  frequencies  could  be  reduced  to  observed  by  addition  of  reasonable 
amounts  of  sediment  to  the  instrument  mass,  as  follows:  a  2-cm-thick 
laver  surrounding  the  roughly  spherical,  buried  NIT ;  a  3-cm-thick  layer 
under  the  flat,  circular  I’L;  two  hemispheres  of  radius  equal  to  that  of  the 
(horizontal)  hearing  area  of  the  buried  spike.  The  reader  is  referred 
to  Appendix  F  tor  an  alternative  investigation  ol  these  parameters. 


a  oia 


s^gma 


S13SI330 


1980  (Appendix  E)  . 
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TABLE  1V-1.  OBS  Coupling  Parameters 


(2) 

(3) 

(4) 

(5) 

(6) 

OBS 

c 

max 

V  / R (s~  ^ ) 
s 

F 

m 

f  (Hz) 
c 

f  (Hz) 

0 

ND 

0 

56 

2.7 

17 

14 

PL 

.53/. 27 

18 

0.22 

19 

14 

SP 

.17 

100 

5.1 

22 

14 

BIO 

.25 

29 

2.2 

9.9 

9.5 

HIGB 

.41 

120 

8.9 

20 

17 

HIGS 

.30 

67 

28 

7.1 

8 

LDGO 

.40 

120 

2  70 

3.6 

12 

MITE 

.53 

50 

1.4 

21 

22 

OSU 

.21 

26 

2.5 

8.1 

10 

sioa) 

.  19 

56 

6.9 

11 

11 

UCSB  l) 

.16 

65/100 

10/36 

10/8.3 

8  1/2 

USGS 

.21 

48 

8.1 

8.4 

8-9 

UTG  1  1 ) 

.23 

200 

100 

9.9 

8  1/2 

UW  1 } 

.  11 

380 

1700 

4.6 

8 

UWF 

.44 

50 

42 

3.9 

8  1/2 

(1) 

(2) 

(3) 

(4) 
(3) 
(6) 


OBS's  having  three  foot  pads; 
Calculated  from  equation  (3) 


VLIICIO  l  la  v  t 


.1/2 


=  10  m/s;  R=(area  of  one  footpad/m)^ 

3  3 

;  q=number  of  foodpads;  p  =1.57  g/cm 


F  =  M/r:  p  RJ; 
m  I  s’ 

f  =  (r  Ru  /4  Mt) 
c  si 


1/2 


V 

s  s 


f  =  observed  coupling  frequency  estimated  from  spectra  of  vertica 
transient  tests 


v  *. 
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The  Q  of  the  coupling  resonance,  theoretically,  is  proportional  to 
Fm  •  Fm,  the  mass  factor,  is  listed  in  Table  IV-1;  a  comparison  with 
the  value  for  Qc  shown  in  Figure  IV-3  indicates  rough  agreement  with 
expectations.  Transient  test  results  in  Appendix  C  can  be  used  for  more 
thorough  comparisons. 

Signals  from  controlled  source  and  transient  tests  can  be  compared 
in  both  the  time  and  the  frequency  domains.  Some  of  the  comparisons  are 
discussed  here  (see  also  Appendix  E)  and  the  remainder  are  presented  in 
Appendix  C.  Data  are  processed  and  presented  in  as  similar  a  format  as 
possible  to  aid  in  comparison.  The  reader  is  referred  to  Appendix  C  for 
the  format  of  the  figures. 

Some  may  question  what  is  learned  from  the  transient  test,  in  that 
it  does  not  force  the  instrument  to  move  in  the  same  way  as  a  seismic 
wave,  and  the  resulting  motions  can  be  complicated.  The  validity  of  using 
the  mechanical  transient  tests  for  obtaining  the  parameters  fc  and  Qc  is 
shown  in  Figures  IV- 5  and  IV-6.  The  time  signals  and  spectra  obtained 
for  the  early  arrivals  from  airgun  shot  182  on  PI.,  SP,  and  SIO  verticals 
are  shown  in  Figure  IV- 5  together  with  the  vertical  transient  tests  of 
SIO.  Note  that  the  spectra  for  PL  and  SP  are  similar,  especially  at 
frequencies  below  20  Hz.  The  SIO  test,  however,  shows  extra  energy  near 
12  Hz  and  a  flattening  of  the  spectrum  above  12  Hz.  The  vertical  transient 
test  for  SIO  (bottom  of  Fig.  IV-5)  shows  a  peak  near  11  Hz;  thus  it  appears 
that  modification  of  the  signal  because  of  coupling  causes  the  spectral 
differences.  A  similar  situation  is  shown  in  Figure  IV-6  for  a  cap  shot. 

Here  the  SIO  test  again  shows  amplification  relative  to  SP  near  11-12  Hz. 

HIGS,  however,  shows  amplification  near  8  Hz.  As  shown  at  the  bottom 
of  Figure  IV-6,  the  HIGS  transient  has  a  peak  near  8  Hz,  again  suggesting 
that  the  effects  of  coupling  are  being  observed  in  the  seismic  response. 

Examples  of  good  coupling  are  shown  in  Figure  IV- 7,  Transient  tests  that 
yield  a  sharp  spike  in  the  time  domain,  such  as  SP  and  MITE,  yield  relatively 
flat  spectral  responses  and  good  seismic  response  to  both  high  and  low 
frequencies.  Two  examples  of  good  coupling  (SP  and  MITE)  and  two  of  poor 
coupling  (USGS  and  UTG)  are  shown  in  Figure  IV-8  for  the  same  CAP  shot. 

(The  low  sensitivity  of  MITE  for  frequencies  below  about  10  Hz  is  of 
internal  origin  and  is  not  related  to  coupling.)  Near  11  Hz  the  USGS 
instrument  shows  a  strong  resonance  and  UTG  shows  a  somewhat  weaker 
resonance.  Spectra  for  both  cap  and  background  noise  for  both  OBS's  fall 
off  above  the  resonance  frequencies.  By  comparison,  the  SP  and  MITE 
verticals  show  relatively  flat  spectra. 

Response  to  Vertical  Input  on  Horizontal  Geophones.  Cross  coupling 
between  vertical  seismic  input  and  horizontal  geophones  does  not  appear 
to  be  a  serious  problem.  If  a  package  is  made  reasonably  symmetrical 
around  a  vertical  axis,  then  horizontal  resonances  will  not  be  excited  by 
vertical  inputs.  Quantitative  data  are  available  from  vertical  transient 
tests  on  OBS's  containing  horizontal  component  sensors,  but  they  have  not 
been  analyzed. 

Horizontal  Input  Recorded  on  Vertical  Geophones.  A  particularly 
good  example  of  this  type  of  cross  coupling  can  be  seen  in  Figure  IV-9. 
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•'ig.  IV-5.  Time  ;ind  frequency  domain  recordings  of  airgun 
shot  182.  The  format  for  these  figures  is  explained  in 
Appendix  C. 
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Fie,.  IV- 6  .  Record!  Hi1,*  shoving  the  correlation  he  tween 
seismic  response  and  transient  tost  response. 
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In  this  1  igure  the  signals  from  all  components  recorded  for  airgun 
shot  182  are  displayed.  The  horizontal  geophones  shown  at  the  bottom 
show  a  strong  6-Hz  arrival  from  about  1.5  to  5  sec  after  the  first 
arrivals.  The  vertical  standards  and  vertical  MITE  and  HIGB  sensors  show 
very  little  of  this  arrival,  indicating  that  it  is  horizontally  polarized. 
Vertical  geophones  from  other  packages,  however,  show  varying  amounts  of 
signal  arriving  during  this  period  which  has  been  transformed  by  the  instru¬ 
ments  from  horizontal  motion  into  vertical  motion.  This  type  of  cross 
coupling  is  undesirable  in  that  it  would  be  extremely  hard  to  predict,  and 
it  represents  a  severe  distortion  of  true  ground  motion.  This  problem  is 
discussed  more  thoroughly  in  Appendix  H. 

The  horizontal  transient  test  recorded  on  the  vertical  geophones 
indicate  the  possible  severity  of  this  problem.  Figure  IV-10  (also 
Figs.  5-10  in  Appendix  C)  shows  results  of  these  tests.  Note  that  in 
the  horizontal  tests  a  second  signal  is  observed.  This  signal  is  believed 
to  be  generated  when  the  electromagnet  armature  iiits  the  pulley  rod 
(Fig.  I  1-5). 

In  Figure  IV-10  the  PI.  standard  shows  a  relatively  sharp  response 
and  flat  spectrum,  indicating  that  cross  coupling  should  not  severely 
affect  this  instrument.  BIO,  however,  shows  a  very  strong  resonance  near 
b  Hz.  This  resonance  indicates  a  possibly  severe  hor izontal-to-vertical 
cross-coupling  problem,  which  is  indeed  seen  in  Figure  1V-9.  The 
transient  tests  for  FCSB  are  interesting  in  that  the  tests  were  done  at 
two  different  levels  rather  than  two  different  azimuths,  and  the  resulting 
transients  excited  two  different  resonances  that  differ  in  frequency  by  a 
factor  of  2  (4  and  8  Hz).  Note,  however,  that  FCSB  did  not  exhibit  a 
severe  cross-coupling  problem  in  Figure  TV-9,  indicating  that  the  wave 
train  observed  in  Figure  IV-9  did  not  contain  sufficient  energy  at  4  or 
8  Hz  to  excite  either  of  these  resonances. 

Response  to  Hor izontal  Input  on  Horizontal  Geophones.  Transient 
tests  were  conducted  in  two  horizontal  directions  or  at  two  different 
elevations  on  each  package.  The  resulting  signals  are  shown  in  Figures  11 
to  13  in  Appendix  C.  In  most  cases  a  resonance  was  observed  that  indi¬ 
cates  motion  in  at  least  one  horizontal  mode,  either  lateral  sliding  or 
rotation  around  a  hor  izontal  axis.  The  complexity  reflects  both  the 
difficult’,'  in  obtaining  reliable  horizontal  transient  data  and  the  true 
complexitv  of  the  motion  of  the  package,  each  with  its  own  peculiarities 
and  symmetries  that  modify  this  type  of  coupling. 


V.  NOISE  MEASUREMENTS 

Two  tvpes  of  noise  studies  were  made:  a  description  of  the  ambient 
background  noise  and  its  mode  of  propagation  and  an  evaluation  of  noise 
induced  in  the  OBS's  hv  water  currents. 

To  describe  the  ambient  background  noise  we  used  the  frequency 
spectra,  the  particle  vector  motion  as  determined  trom  the  three  components 
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.  I V  — ! 0 .  Horizontal  transient  tests  recorded  on  vertical 
geophones .  The  severity  of  this  type  of  cross  coupling 
foi-  some  instruments  is  obvious.  From  Sutton  et  al.,  1980 
(Append  i:.  K)  . 


of  motion,  and  the  frequency-wave  number  spectra  determined  from  the  array 
of  DBS's.  We  found  that  the  noise  spectra  peaked  very  strongly  around 
2.5  to  3  Hz  and  in  this  frequency  band  the  particle  motion  was  elliptical 
and  had  phase  velocities  of  about  20  to  50  m/sec;  these  parameters  are 
consistent  witli  the  propagation  of  the  noise  as  surface  waves.  We  also 
found  that  shot-generated, dispersed  surface  waves  have  properties  similar 
to  the  background  noise,  supporting  the  contention  that  the  mode  of 
generation,  propagation,  and  attenuation  of  the  noise  can  be  understood 
in  terms  of  dispersed  surface  waves.  Further,  and  possibly  of  greatest 
practical  importance,  we  established  that  these  surface  waves  are  highly 
attenuated  in  the  water  and  weakly  recorded  on  a  hydrophone.  This  implies 
that  hydrophones  in  deep-sea  sediment  environments  are  likely  to  be 
substantially  quieter  and  give  less  complicated  seismograms  than  geophones. 
The  determination  of  current- induced  noise  on  the  OBS's  was  evaluated  by 
comparing  the  spectra  of  several  instruments  against  current  speed.  The 
results  from  this  comparison  are  not  definitive. 

A  more  detailed  discussion  of  these  results  is  given  below. 


Character ist ics  of  Ambient  Seismic  Noise 

An  example  of  ambient  noise  and  its  amplitude  spectrum  as  recorded 
on  the  plate  vertical  is  shown  in  Figure  V-l,  i rom  record  210.  The 
spectrum  shows  that  the  noise  is  dominated  by  energy  in  the  2-  to  A~Hz 
band.  A  vector  plot  of  the  particle  motion  is  shown  in  Figure  V-2.  In 
these  plots  the  vertical  axis  is  one  component  of  motion  and  the  hori¬ 
zontal  axis  is  an  orthogonal  component.  In  these  graphs  rectilinear,  or 
body-wave  motion,  would  appear  in  the  vertical-horizontal  plots  as  vectors 
of  constant  slope,  whereas  elliptical  particle  motion  would  appear  as 
fan-shaped  patterns.  It  is  clear  from  Figure  V-2  that  the  predominant 
motion  is  elliptical.  Another  example  is  noise  preceding  airgun  shot  182. 
The  three  components  of  motion  as  recorded  on  the  plate  and  the  associated 
particle  motion  are  shown  in  Figures  V-3  and  A. 

rn  this  experiment  the  two-dimensional  array  can  he  used  to  determine 
the  coherency  and  phase  velocity  of  the  noise  by  taking  the  three- 
dimensional  Fourier  transform,  or  frequency-wave  number  spectrum  of  the 
arrav  data.  An  example  of  the  wave  number  spectrum  taken  at  2  Hz  is 
shown  in  Figure  V-3.  From  this  spectrum  the  phase  velocity  at  2  Hz  is 
found  to  he  about  30  m/sec.  For  a  more  detailed  discussion  of  tills 
analysis  see  Appendix  0. . 


Comparison  of  Noise  witli  S ho L-fA'inyratcyd  Sur ta_ce_W  ayes 

Nearly  all  the  cap  and  airgun  shots  generated  low-velocity  surface 
waves  in  the  2-  to  A-Hz  band.  An  example,  airgun  shot  j-82,  as  recorded 
bv  tiie  arrav  is  shown  in  Figure  V-h.  (See  also  iig.  IV-9.)  lhe  spectrum 
ol  this  wave  as  recorded  on  the  plate  vertical  for  airgun  shot  131  is 
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Noise  record  preceding  airgun  shot  182  showing  vertical  and  horizontal  components  of  motion 
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shown  in  Figure  V-7  and  is  very  similar  to  the  noise  spectrum  in  Figure  V-l. 
i'lie  particle  motion  tor  this  wave  packet  is  elliptical,  as  shown  in 
Figure  V-8 .  A  group  velocity  and  phase  velocity  dispersion  curve  for  this 
wave  packet  is  shown  in  Figure  V-9  (from  Tuthill  et  al.,  1980  ,  Appendix  G) . 

The  similarity  of  the  spectra,  dispersion  curves,  and  particle 
mot  ion  lor  the  noise  and  the  shot-generated  surface  waves  strongly 
suggests  that  the  noise  is  propagating  as  short-wavelength  dispersive 
surface  waves.  The  wavelengths  in  the  2-  to  4— H z  band  are  found  from  the 
corresponding  phase  velocities  to  he  about  11  m  to  6  m. 


Comparison  Between  Hydrophone  and  Se  i  smome  ter  j)a_t  a 

Figure  V— 10  shows  a  comparison  between  the  LUGO  hydrophone,  the 
plate  vertical  and  the  plate  horizontal  for  airgun  shot  182.  It  is  clear 
from  tliis  figure  that  the  surface  wave  train  is  only  weakly  recorded  on 
the  hydrophone.  This  is  to  he  expected  since  this  motion  is  caused  by 
short-wavelength  surface  waves  that  are  highly  attenuated  in  the  water. 

The  previous  discussion  also  implies  that  the  2-  to  4-Hz  noise  should  be 
greatly  reduced  on  the  hydrophone.  This  implication  is  verified  in 
Figure  V-ll,  a  comparison  of  noise  spectra  from  the  LOGO  hydrophone  and 
the  IlIGS  vertical.  This  result  alone  implies  that  hydrophone  records 
should  he  quieter  and  simpler  than  geophone  records  because  they  will  be 
contaminated  bv  short-wavelength  ambient  noise  and  short-wavelength  signal- 
i minced  noise.  This  is  particularly  important  in  crustal  refraction 
protiles  because  wavelengths  of  about  10  m  will  contain  information  about 
tiie  immediate  neighborhood  of  the  instrument  only  and  is  of  no  use  in 
determining  average  crustal  properties.  Worse  than  this,  the  short- 
wavelength  noise  will  greatly  degrade  the  long-wavelength  signals. 


Correlation  of  Background  Noise  w it h  Currents 

Tlu’  purpose  of  this  section  is  to  investigate  the  possibility  that 
(a)  currents  cause,  shaking  or  the  instruments  and  hence  noise;  or  alterna¬ 
tive  1  v  ,  (h)  that  water  currents  traversing  the  array  couple  noise  into 
the  bottom  tiie  propagates  across  the  array. 

During  one  13-hour  period,  hourly  samples  (oi  approximately  1-min. 
duration)  of  background  noise  were  digitized.  Water  currents  recorded 
on  each  of  the  four  current  meters  for  15  min  preceding  and  following 
each  noise  sample  were  averaged  and  ranked  and  are  plotted  in  Figure  V-12. 
During  this  time  period  the  currents  ranged  I rom  1  to  4-1/2  cm/sec.  This 
range  is  about  as  great  as  was  observed  during  the  entire  experiment.  The 
larger  standard  deviations  for  the  higher  current  levels  suggest  a 
greater  amount  of  irregularity  in  current  across  the  CBS  array  during 
those  times. 

Noise  spectral  estimates  were  obtained  for  each  1-min  sample  by 
dividing  the  record  into  N  4— sec— long  records  (whore  N1  30)  and  averaging 
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-10.  Comparison  of  recordings  by  the  plate  and  the  LOGO  hydrophone  of  airgun  shot  182.  Note  the 
h  attenuation  of  the  surface  waves  and  horizontal  motion  on  the  hydrophone. 
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the  N  spectra.  For  each  frequency  the  mean  X  and  the  2-sigma  confidence 
limits  were  obtained.  Tn  addition,  the  normalized  variance,  defined  by 


(Xi  -  X)2 


i  =  1 


was  determined.  For  a  stochastic  or  random  process  the  variance  will  be 
proportional  to  the  mean,  who r as  for  a  deterministic  process  the  variance 
will  he  zero.  An  additional  statistic,  the  fractional  standard  deviation, 
was  also  calculated  and  is  dol'incd  by 


1  N  . 
l  i  = 


(Xi  -  X) 


For  a  gauss ian  process  this  statistic  will  bo  equal  to  a  constant,  whereas 

for  a  deterministic  process  this  statistic  will  tend  to  zero. 

Examples  of  these  spectra  and  their  statistics  are  shown  in  Figure  V-13 
1  or  the  Spike  V  and  the  I’CSB  V.  It  can  be  seen  from  these  figures  that 

Luo  hO-Hz  noise  and  noise  peaks  between  20  and  10  Hz  have  low  standard 

errors,  indicating  that  they  are  deterministic,  as  expected  for  the  60-Hz 
noise.  Also,  the  20  to  30  Hz  may  have  a  mechanical  source  (for  example, 
i rom  a  motorboat.)  The  predominant  noise  at  2  Hz  is  more  nearly 
stochastic. 


interpretation  of  these  spectra  is  made  difficult  because  of  the 
several  possible  ways  in  whirl)  currents  could  induce  noise.  For  example, 
currents  might  shake  one  instrument  (such  as  the  moored  current  meters) 
more  than  the  others  and  this  energy  could  be  coupled  to  the  bottom  and 
radiate  across  the  array  as  surface  waves  in  t ho  2-  to  4-Hz  band,  or  the 
shaking  of  individual  instruments  might  stimulate  coupling  resonances 
observed  in  the  lift  tests  and  these  resonances  would  occur  at  frequen¬ 
cies  peculiar  to  each  OHS. 

io  test  whether  current-related  noise  is  coupled  into  the  ground 
and  radiated  across  the  array  we  have  plotted  the  height  of  the  2-  to 
4-!lz  spectral  peak  on  the  Spike  V  versus  current  speed  in  Figure  V-14. 
the  correlation  between  these  two  parameters  is  low,  suggesting  that 
there  is  Little  cur  rent- i nduced  noise  propagating  across  the  array.  (Note 
that  the  Spike  V  was  chosen  for  this  test  because  it  was  buried  and  should 
in  insensitive  to  current- induced  slinking.) 

I’d  test  whether  individual  instruments  were  sensitive  to  current 
slink  in-,  we  used  as  an  example  tin'  IH'.SB  V  data. 


average  of  13  spectra 
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It  can  be  seen  on  Figure  V—  Hi  that  the  UCSBV  lias  a  noise  peak  at 
about  9  Hz  that  is  not  seen  by  the  Spike  V.  The  horizontal  pull  lift 
tests  showed  this  instrument  to  have  a  strong  resonance  at  this  frequency 
and  it  is  clear  that  this  resonance  is  being  induced  by  the  noise.  To 
evaluate  whether  this  instrument  was  sensitive  to  the  currents  we  plotted 
the  amplitude  of  the  2-  to  4-Hz  and  9-Hz  noise  peaks  versus  current  speed 
in  Figure  V-19.  The  lack  of  correlation  suggests  that  the  currents  are 
not  causing  the  noise. 

We  conclude  that  in  this  range  of  current  speed  11  to  4  cm/sec)  the 
currents  are  not  the  source  ot  the  noise.  At  this  time  we  have  no  idea 
what  is  causing  the  noise. 

Additional  noise  data  are  presented  in  Appendix  I). 


VI.  CONCLUSIONS 


On  the  basis  of  the  preceding  discussions  we  arrive  at  recommendations 
lor  the  design  criteria  for  ocean  bottom  seismic  systems.  Although  none 
oi  these  recommendat ions  is  absoluLel'  necessary  for  obtaining  high- 
I  idelitv  seismic  signals  at  any  particular  site,  if  the  recommendations 
arc  followed,  OBS  svstems  can  he  built  that  will  be  able  to  obtain  higli- 
lidelitc  seismic  signals  anywhere. 

Tiie-  primary  purpose  of  an  OBS  is  to  record  the  motion  of  the  ground. 

If  a  record  of  this  motion  cannot  be  obtained  with  reasonable  fidelity 
because  oi  cither  (1)  high  noise,  (2)  poor  coupling,  (3)  lack  of  dynamic 
range,  (4)  lack  of  storage  capaeitv,  or  (5)  lack  of  reliability,  then  the 
DBS  lias  at  least  partially  failed.  This  test  addressed  the  first  two  of 
these  problems.  Indeed,  it  coupling  and  noise  problems  are  serious,  no 
amount  of  dynamic  range  or  storage  capability  will  help. 

Several  observations  from  the  Lopez  Experiment  are  important  when 
considering  mod i t ications  to  present  systems  or  design  of  new  OBS's; 

(1)  Tlie  geophone  standards  produced  generally  consistent  results  and  seemed 
to  record  true  ground  motion.  A  design  goal  might  he  to  build  an  OBS 
similar  to  one  of  the  standards — either  plate-like,  with  neutral  density 
in  the  sediment,  or  shoved  into  the  sediment.  None  of  the  standards 
however ,  was  built  to  record  data  internally,  free-fall  to  the  bottom,  and 
return  to  the  surtace.  Adding  hydrodynamic  stability  during  drop, 
buovanev  tor  return,  and  recording  capability  would  change  the  basic 
features  ot  these  units  unless  the  sensors  wore  caret ully  located  in  the 
package  or  median  i  ca  1 1  v  decoupled  I  rotn  the  record  i  tg— retrieval  unit. 

Indeed,  the  M LIT  and  I1IDB  instruments,  which  showed  good  coupling  character- 
j.,t  ics,  have  c ha rac t or  is t  ics  similar  to  the  plate  standard.  An  OBS  with 
character ist  irs  similar  to  the  neutral  density  or  spike  standards  would  not 
be  suitable  for  hard-bottom  applications,  nor  world  emplacement  and 
retrieval  be  a  simple  matter  in  soft  sediment. 


30-000 


2-000  4.000  6-000  8-000  10-000  12-000 

CURRENT  INDEX 


NOISE  VS  CURRENT 


2-000  4-000  6-000  8-000  10-000  12-000 

CURRENT  INDEX 

NOISE  VS  CURRENT 


.  Plot  O  i  the  .imp!  i  t tide  of  the  '3-Hz  and  9-11?.  noise  peaks 

v.-r.  us  i  nrretit  index  for  the  PCS  I!  V.  There  is  not  a  strong  eor- 
r--l.it  imi  between  enrrent  speed  and  noise  levels.  The  current 
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(2)  Hydrophones  yield  records  that  are  not  distorted  by  coupling 
problems  or  complicated  by  low-velocity  shear  arrivals.  However,  the 
tailure  to  record  (near  normally  incident)  shear  arrivals  can  be  a  curse 
as  well  as  a  blessing  in  some  applications.  This  is  easily  seen  in 
Figure  IV-9,  where  the  hydrophone  sees  very  little  signal  after  the  first 
arrivals.  Although  hydrophones  may  be  adequate  for  most  refraction  work 
and  are  useful  for  earthquake  studies,  the  loss  of  shear  wave  data  in 
either  application  can  be  considered  a  serious  shortcoming.  Hydrophones 
should  he  considered  a  valuable  sensor  in  (IBS's,  but  not  a  replacement 
for  geophones. 

(1)  Maximizing  base  dimensions  relative  to  height  of  the  package 
appears  to  reduce  cross  coupling  caused  by  rocking.  Placing  geophones 
near  the  center  of  cross-coupling  rotation  should  also  lessen  this  effect, 
especially  for  the  vertical  component.  Reducing  the  cross  section  in 
water  would  certainly  also  lessen  possible  problems  with  bottom  currents. 
Although  current-generated  noise  was  not  observed  at  Lopez,  the  highest 
currents  observed  during  the  experiment  are  often  exceeded  in  the  deep 
ocean  (Kasahara  et  al.,  1979).  There  is  some  published  evidence  for 
correlation  of  current-generated  noise  with  OBS  package  configuration 
(Bradner  et  al.,  1969  •,  Kasahara  et  al.,  1979  )  Kasahara  et  al.,  in  press  ; 
and  Duennebier  et  al.  ,1980.)  Intuitively  one  would  expect  a  squat, 
low-profile  instrument  design  to  he  preferable;  this  intuitive  argument 
may  extend  even  to  the  flotation  and  recording  packages  of  the  instruments 
with  externally  deployed  sensor  packages. 

(4)  First  order  theory,  assuming  a  damped  harmonic  oscillator, 
suggests  that  increasing  the  surface  area  in  contact  with  the  sediment 
increases  the  coupling  frequency  and  the  damping.  This  effect 
forces  the  DBS  motion  to  be  closer  to  that  of  the  sediment.  In  addition, 
more  sediment  and  water  are  likely  to  be  entrained  with  the  OBS,  thus 
causing  a  drop  in  the  coupling  coefficient,  C,  and  improved  response  to 
frequencies  above  the  coupling  frequency.  In  general,  the  Lopez  observa¬ 
tions  support  tli  is  theorv. 

The  incorporation  of  an  in  situ  transient  test  device  on  existing 
(IBS's  appears  feasible.  In  addition  to  providing  a  good  check  on  overall 
system  operation,  it  would  give  two  of  the  three  principal  coupling  para¬ 
meters,  tc  and  Qc,  needed  for  calculating  coupling  response.  The  third, 
the  coupling  coefficient,  C,  can  be  estimated,  as  discussed  in  section  IV. 

(9)  Adjusting  the  density  of  the  package  to  be  slightly  greater  than 
that  of  water  also  would  decrease  the  value  of  C,  decreasing  the  effects 
of  coupling.  This  would  make  the  instrument  more  sensitive  to  high 
frequencies  and  reduce  the  resonant  amplification.  Special  care  must 
be  taken  in  this  case  to  keep  a  low  profile  and  good  contact  with  the 
bottom  to  reduce  cross  coupling  and  current  problems. 

One  problem  not  discussed  in  this  report  is  horizontal-to-horizontal 
coupling.  Unlike  vertical-to-vert ical  coupling,  it  is  much  more 
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susceptible  to  torques  and  asymmetries  In  package  configuration  and 
emplacement.  Asymmetry  can  cause  considerable  azimuthal  variation  in 
response  to  horizontal  motion.  The  complexity  of  the  horizontal  coupling 
problem  makes  it  difficult  to  understand,  although  following  basic  princi¬ 
ples  such  as  keeping  asymmetries  to  a  minimum  and  decreasing  the  cross 
section  in  the  water  should  help  considerably. 

In  the  l.opez  tests  the  amplitude  responses  of  the  spike  and  plate 
standard  instruments  were  generally  within  2  dB  of  each  other  and  the 
waveforms  were  well  duplicated.  We  judge  from  the  similarity  of  the 
responses  that  these  instruments  were  responding  accurately  to  the  motion 
of  the  upper  few  centimeters  of  the  sediment  layer.  We  believe  that  we 
now  know  enougli  about  problems  with  noise  and  coupling  to  design  an  OBS 
system  that  will  sense  true  ground  motion  with  fidelity  equal  to  that  of 
the  standards.  We  also  know  that  in  some  places  geophones  will  be 
particularly  noisv  because  of  short-wave  length  signals  in  the  sediments, 
in  those  areas,  particularly,  the  addition  of  a  hydrophone  to  the  OBS  is 
highly  desirable. 

Knowing  how  to  sense  the  motion  is  only  one  problem.  Recording  and 
retrieving,  the  data  deserve  equal  study  and  effort.  Numerous  technological 
opt  ions  are  available:  analog  or  digital  recording;  acoustic  or  timed 
releases;  spherical  or  cylindrical  pressure  cases;  motor  or  explosive  bolt 
releases  and  so  forth.  I'nl  ike  Lite  problems  with  coupling  and  noise,  we 
ran  anticipate  that  more  and  bettor  solutions  will  become  available  for 
these  problems  as  technology  advances. 


VI  1  .  Kb'TIcRK  DEVELOPMENTS 


The  coupling  and  noise  problems  in  existing  instruments  can  be 
reduced  significantly  by  modifications  suggested  by  the  results  of  this 
experiment.  With  such  modifications  these  instruments  will  be  adequate 
to  continue  fruitful  research  in  marine  seismology  over  the  next  few  years. 
This  research  not  only  will  result  in  significant  progress  in  seismology, 
hut  also  will  provide  additional  valuable  data  on  design  parameters, 
particularly  if  collaborative  field  work  continues  in  which  the  performance 
of  instruments  oi  different  design  can  he  compared. 

The  above  notwithstanding,  new  instrumentation  should  be  designed 
and  tested.  Ocean  bottom  seismology  is  far  too  new  a  field  to  permit 
stagnation  in  the  development  of  instrumentation. 

three  systems  appear  desirable  for  the  foreseeable  future:  (1)  a 
simple,  inexpensive  ^hydrophone  '? )  system  that  will  allow  use  of  a  large 
number  of  units  for  dense  arrays;  (2)  a  somewhat  more  expensive  system 
with  a  hydrophone  and  three  geophones  to  augment  arrays  ot  the  simpler 
systems,  to  he  used  in  detailed  micro-earthquake  studies,  and  (3)  a 
relatively  expensive  system  capable  of  handling  a  variety  of  sensor  systems 
including  hydrophone  arravs,  geophones,  and  long-period  sensors. 

Tlie  simple  arrays  could  be  easily  deployed  and  retrieved  and  would 
be  cheap  enough  so  that  large  number  of  units  could  be  used  in  any 
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expor  imi-fiL .  Those  units  would  ho  of  primary  valuo  in  refraction 
experiments  and  couid  he  used  to  provide  dense-  coverage  in  micro-earthquake 
s  t  ud i es . 


file  tour  component  units  would  supply  sediment  and  shear  velocity 
information  for  refraction  lines  and  shear  wave  arrival  times,  particle 
motion,  and  surtace  wave  information  in  earthquake  studies.  These  units 
would  require  more  care  in  emplacement  and  recovery  and  would  be  more 
costly.  The  "expensive  system"  would  he  t  lie  experimental  unit  used  in 
special  studies  and  for  test  ini;  o|  new  ideas. 

Certain  ocean  bottom  se  i  smo  1  og  i  ca  1  research  is  also  expected  to 
require  instnumnts  capable  ol  vet-.  Ion.;  deployment  times.  Development 
in  this  area  also  should  he  supported. 
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CATALOG  OF  LOPEZ  ISLAND  EXPERIMENT  DIGITIZED  EVENTS 


On  the  following  pages  is  a  catalog  of  the  seven  digital  tapes 
produced  by  the  SIO  Geospace  recording  system.  Each  component  of  each 
instrument  is  identified  by  a  code  number  described  on  the  next  page. 
Catalog  entries  are  listed  according  to  Geospace  record  number.  The 
code  numbers  for  all  components  recorded  during  an  event  are  listed 
next  to  each  Geospace  record  number. 

Record  number  500  was  not  recorded  digitally  and  was  hand-digitize 
from  a  paper  recording.  Records  146  to  148  were  not  recorded  correctly 
Data  for  these  events  were  re-digitized  from  the  DIGS  analog  tapes 
recorded  inside  the  OBS.  These  records  are  not  in  this  index  but  are 
found  elsewhere  in  the  report  (Appendix  E,  Fig.  E-5) . 
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LIST  OF  LOPEZ  COMPONENTS 

SIO 
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Number 

Digital 
Log  Book 

Present 

Notation 

Descript  ion 

1 

SP  1 

SP  V 

Spike 

standard  vertical 

2 

SP  2 

SP  N 

Spike 

standard  horizontal 

3 

SP  3 

SP  E 

Spike 

standard  horizontal 

4 

ND  1 

ND  V 

Neutral  density  standard  vertical 

5 

ND  2 

ND  N 

Neutral  density  standard  horizontal 

6 

ND  3 

ND  E 

Neutral  density  standard  horizontal 

7 

PL  1 

PL  V 

Plate 

standard  vertical 
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PL  N 

Plate 

standard  horizontal 
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PL  E 

Plat  e 

standard  horizontal 
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Un  i  v  . 

of  Washington  vertical 

11 
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Uni  v. 

of  Washington  horizontal 

12 
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U  n  i  v  . 

of  Washington  Flower  Pot  vert 
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Scripps  Institute  vertical 

14 
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Bedford  Institute  vertical 
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Bedford  Institute  horizontal 
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State  vertical 
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SELECTED  DIGITIZED  EVENTS  RECORDED 
DURING  THE  LOPEZ  ISLAND  EXPERIMENT 


SELECTED  DIGITIZED  EVENTS  RECORDED  DURING  THE  LOPEZ  ISLAND  EXPERIMENT 


This  appendix  contains  111  records  of  signals  recorded  on  ocean  bottom 
seismometers  at  Lopez  Island.  Each  recording  was  processed  identically 
within  the  bounds  discussed  below.  Records  are  presented  four  per  page 
with  time  domain  recordings  on  the  left  and  frequency  domain  on  the  right. 
Each  time  domain  recording  is  1024  samples  long;  the  iiorizontal  line  under 
each  recording  is  2  sec  (500  samples)  long.  The  vertical  lines  under  each 
record  show  the  section  of  the  recording  that  was  transformed  into  the 
frequency  domain  (power  spectrum)  shown  on  the  right.  In  each  case,  the 
section  to  be  transformed  was  convolved  with  a  104  cosine  taper  window 
and  placed  in  the  middle  of  a  512-point  window  after  decimation  by  a  factor 
of  2.  No  corrections  are  made  for  possible  aliasing.  The  resulting  512- 
point  time  series  was  then  fourier  transformed  and  the  "periodogram"  was 
computed  to  obtain  an  estimate  of  the  power  spectrum.  In  addition  to  the 
signal,  a  section  of  noise  (the  same  length  as  the  signal)  was  transformed 
and  plotted  (lighter  line).  This  makes  it  possible  to  distinguish  energy 
generated  b>  signal  from  noise.  The  vertical  scale  for  each  spectrum  is 
dH  referenced  to  1  digital  unit.  The  horizontal  scale  is  in  Hertz.  Under 
each  spectrum  is  the  information  necessary  to  identify  each  event.  The 
first  two  to  four  letters  are  the  instrument  identifier  used  throughout 
this  report  followed  by  the  component  designation  (V,  H,  or  P,  for 
vertical,  geophone,  horizontal  geophono,  or  pressure).  Then  given  are  the 
type  of  source.  Geospace  record  number,  and  channel  number  from  Appendix  B. 
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NOISE  STUDIES 


The*  figures  in  this  appendix  are  discussed  in  section  V  of  this 
report.  Figures  D-l  to  !>5  were  computed  from  the  hourly  noise  samples  of 
dune  25,  l'7S.  The  spectra  are  listed  according  to  the  current  rank  from 
Figure  1 2 .  The  number  1  represents  a  noise  sample  from  the  hour  in 
which  the  weakest  current  was  measured  and  12  corresponds  to  the  strong¬ 
est.  errent.  Tin  spec’.  ra  «.re-  labeled  with  the  sensor  name  and  digital 
record  :.-unKei  '  rom  Appendix  B  For  1.DC0  and  H1GS  components,  noise 
simples  were  -.ttken  hiring  the  last  four  sampling  periods  only. 

Details  of  the  ipoctral  an  r lysis  »ie  given  here.  The  hourly  noise 
i-afip  i  .ere  ,ib  -ut  .  or  2  min  long.  From  these  samples,  the  "most  quiet" 
id.'1’-  Jh  ii  .  i  -..ample  blcck  (about  4  sec)  was  decimated  to  512  sample  points 
and  win...'";!  i-ich  j  1 02  cosine  window.  A  Fast  Fourier  Transform  was  used 
to  prod’i  e  a  2 5b -point  power  spectrum  with  a  Nvquist  frequency  of  62.5  Hz. 
Ih'*  pov-o  spectral  estimates  are  plotted  in  dB  referenced  to  one  digital 
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G.  H.  SutCon,  F.  K.  Duennebier,  and  B.  Iwatake 

Hawaii  Institute  of  Geophysics 
University  of  Hawaii 
Honolulu,  Hawaii  96822 


Abstract 


Unlike  response  of  seismometers  resting  on  hard  rock  where  the 
seismometer  case  moves  with  the  rock  to  high  frequencies,  the  response  of 
ocean  bottom  seismometers  (OBS)  can  be  strongly  affected  by  the  low  me¬ 
chanical  strength  of  ocean  sediments.  The  motion  as  measured  bv  the 
seismometer  will  not  ft  ow  the  expected  relationships  between  pressure 
and  particle  motion  for  different  wave  types.  Cross  coupling  between 
horizontal  and  vertical  motions  can  occur,  especially  when  there  is 
differential  motion  between  water  and  sediment.  Resonant  amplification 
and  attenuation  of  higher  frequencies  also  occur.  Secondary  seismic 
arrivals  are  especially  subject  to  distortion.  Overall  response  is  strongly 
dependent  upon  the  mass  and  configuration  of  the  OBS  and  the  rigidity  and 
density  of  the  bottom  material.  Tests  at  Lopez  Island,  Puget  Sound  using 
both  directly  applied  mechanical  transients  and  seismic  signals  with 
various  instrument  configurations  demonstrate  the  above  effects  and  provide 
some  guidance  for  improved  designs. 
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Introduction 


The  usual  reason  for  placing  a  seismometer  anywhere  is  to  measure 
free-field  particle  motion  in  response  to  some  seismic  source.  The 
presence  of  the  seismometer  does  not  usually  modify  the  particle  motion  at 
seismic  frequencies  if  care  is  taken  to  emplace  the  instrument  on  a  rela¬ 
tively  rigid  spot.  On  the  ocean  floor,  however,  the  most  common  plant  is 
thick,  soft,  low-velocity  sediment.  For  this  case,  the  presence  of  the 
seismometer  can  severely  distort  the  free-field  motion.  In  addition,  the 
seismometer  is  emplaced  at  a  boundary  between  a  liquid  and  a  solid,  unlike 
land  emplacement,  and  the  effects  of  buoyancy  can  strongly  affect  the 
observed  motion.  The  boundary  also  complicates  the  situation  because 
strong  differential  motion  can  occur  across  it. 

The  problem  to  determine  the  response  of  an  OBS  to  bottom  motion  is 
similar  to  that  faced  a  number  of  years  ago  by  geophysicists  concerned 
with  the  response  of  prospecting  geophones  placed  on  compliant  materials, 
with  the  additional  complication  of  water  motion  and  buoyancy  forces 
(Washburn  and  Wiley,  1941;  Wolf,  1944;  Bycroft,  1956;  Lamer,  1969). 

For  the  prospecting  case,  the  response  depends  principally  on 
rigidity  and  density  of  soil  and  mass  and  bearing  area  of  geophone.  Damp¬ 
ing  results  from  reradiation  of  seismic  energy  as  well  as  from  internal 
dissipation  (the  latter  is  generally  ignored)  and  depends  upon  the  ratio 
of  instrument  radius  to  a  (shear)  wavelength.  With  low  rigidities  of 
ocean  bottom  sediments,  these  wavelengths  can  be  shorter  than  we  might 
first  guess,  e.g.,  for  V  =  10  m/s  and  f  =  100  Hz,  \  ~  10  cm.  The  curves 
in  Figure  E-l  show  corner  (resonance)  frequency  and  resonance  amplifica¬ 
tion,  Q,  as  functions  of  dimensionless  instrument  mass  and  ratio  of  shear 
velocity  to  instrument  radius.  For  a  given  uniform  bottom  sediment,  corner 
(resonance)  frequency  increases  approximately  as  (R/M)^'^  and  Q  increases 
as  (M/R^)l'“  where  R  is  the  radius  of  the  cross  section  in  contact  with 
the  sediment  and  M  is  the  instrument  mass. 

In  this  paper  we  provide  a  simple  extension  of  the  theory  to  include 
buoyancy  forces  (for  vertical  motion  only)  ,  present  the  theory  for  a 
mechanical  transient  test  designed  to  excite  coupling  resonances,  and 
compare  theoretical  expectations  with  observational  results  obtained  from 
special  tests  in  Shoal  Bay,  Lopez  Island,  Puget  Sound.  Sutton  et  al. 
(1980),  Zelikovitz  and  Prothero  (1980),  Tuthill  et  al.  (1980),  Lewis  and 
Tuthill  (1980),  and  Johnson  and  McAlister  (1980)  present  additional 
results  related  to  bottom  coupling,  based  mainly  on  the  Lopez  Island  OBS 
Intercomparison  Experiment. 


Observations 


Figures  E— 2  and  E-3  show  a  P  arrival  and  spectra  for  arrivals  from  a 
local  earthquake  recorded  in  deep  water  off  the  southeast  coast  of  the 
island  of  Hawaii  on  three  standard  HIG  OBS's  (HIGS,  Sutton  et  al.,  1977a). 


Corner  frequency,  f0  (Hz) 
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Fig.  F.-l.  Coupling  frequency  fQ  and  resonancy  amplification,  Q,  of  a  mass, 
M,  on  the  free  surface  of  an  elastic  half  space.  R  =  radius  of  bearing 
area;  Vs  and  Ps  are  shear  velocity  and  density  in  the  half  space,  respec 
tivelv . 
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The  resonance  peaks  shown  (3.7,  5.3,  and  6  Hz)  are  found  on  all  arrivals 
on  the  respective  OBS's  during  the  same  deployment.  Interpreted  as  a 
coupling  resonance,  the  Q  is  about  7  to  20  and  the  sediment  shear  velocity 
Vs,  is  between  5  and  10  m/s  (see  Fig.  E-l). 

Figure  H-4  is  a  comparison  of  HICS  and  HIGB  recordings  of  the  same 
shot.  The  HIGB  was  developed  to  improve  bottom  coupling  and  reduce  noise 
from  interaction  with  bottom  currents  and  from  the  tape  recorder  (Byrne  et 
al.,  H77).  Tiie  geophones  in  the  HIGB  are  placed  in  a  small  independent 
pressure  case  and  separated  from  the  main  package  by  about  1  m 
(Sutton  et  al.,  1980).  The  OBS's  were  located  within  tens  of  meters 

of  each  other  in  deep  water.  The  hydrophones  are  well  correlated,  showing 
the  similarities  in  wave  forms  expected  for  instruments  next  to  each  other 
There  should  be  no  difference  between  the  hydrophones;  however,  the 
modified  instrument  hvdrophone  was  recorded  at  12-dB  lower  gain  because  of 
a  cross-taLk  problem  with  the  time  code  channel.  All  channels  have  been 
identically  filtered  to  reduce  high-frequency  noise.  The  vertical  geo¬ 
phones  show  the  effects  of  deploving  the  geophones  away  from  the  main 
instrument  package .  The  standard  instrument  (HITS)  gain  was  18  dB  lower 
than  the  modified  (HIGB)  for  both  the  vertical  and  the  horizontal;  thus, 
since  the  gains  are  set  automatically  inside  each  OBS,  the  noise  level  is 
about  8  times  lower  for  the  modified  instrument  than  for  the  standard. 

Note  that  the  frequency  content  and  wave  forms  received  by  the  modified 
instrument  vertical  component  are  similar  to  those  of  the  hydrophone, 
whereas  the  wave  form  of  the  standard  vertical  shows  a  strong  resonance  at 
low  frequency  and  little  coherence  with  the  wave  form  on  the  hydrophone. 

A  teleseismic  event  that  was  clearly  recorded  on  HIGS  and  HIGB 
instruments  during  the  Lopez  experiment  (Sutton  et  al..  1980; 
is  shown  in  Figures  F-5  and  F.-6.  The  instruments  were  3-1/3  m  apart.  At 
270-km  range,  primary  arrivals  should  arrive  within  about  20°  of  vertical 
and  be  of  compress iona L  type  for  the  time-interval  shown.  In  the  early 
part  of  the  record,  the  hydrophone  and  HIGB  (BOOBS)  vertical  show 
excellent  correlation  as  expected  for  near-vertical  incidence.  In  the 
latter  part,  correlation  is  not  as  good,  with  the  vertical  showing  some 
excess  low-frequency  energy  and  similarity  with  the  horizontal  component. 
This  may  result  from  converted  "granddaddy"  type  waves  that  have  a 
relatively  weak  pressure  component  (Tuthill  et  al, ,  1980).  These 
converted  waves  produce  spectral  peaks  near  3  Hz  on  all  the 
components.  The  HIGS  (FOBS)  vertical  shows  an  emergent  character 
an  excess  of  low  frequency  resulting  from  some  combination  of  coupling 
resonances,  high-frequency  attenuation,  and  sensitivity  to  cross-coupling 
from  horizontal  input.  Because  of  the  stiffness  of  the  Lopez  sediments 
the  HIGS  (and  all  the  other  OBS's)  probably  rested  higher  in  the  water 
than  for  normal  deep  ocean  sediment.  The  spectrum  from  the  HIGS  horizon¬ 
tal  shows  a  resonance  near  8  Hz  and  resulting  drop  at  higher  frequencies 
in  comparison  with  the  HIGB  horizontal,  as  expected  from  the  coupling 
theory  discussed  in  the  following  section. 
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Fig.  K-  5 .  Earthquake  (nr  large  quarry  M;W',  M  =  '3.4,  H  =  3.1  km,  A=  270  km)  recorded 
during  Lope’  .  xperirn  ni  on  ilIC.S  ,'.NI)  H It "li  in  water  depth  of  about  8  m. 
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Spectra  ol  signal  shown  in  Figure  1 
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We  derive  nere  the  t  henry  of  coupling  of  vertical  motion  of  the  ocean 
bottom  to  vertical  motion  of  an  CBS-  The  DBS  shape  is  simplified  to  that 
of  a  vertical  cvl  hide  ■'  to  enable  us  to  ignore  complications  from  shape  and 
change  in  bearing  area  v  •  th  buoyancy  and  rigidity  of  the  bottom. 

Keferr  in.;  to  fine.  :  w  ■  del  ir.e  Lae  following  parameters: 

nasi  Lien  e.  j.'.u  .  'l  ee 
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:  ft  .  t  ; > ■  :  ustn-.t  o:  the  sediment-OBS  contact  area 
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=  .‘.it  f  ”M  r.t  :  ■  it  iea  of  idBS 

.  -Ah  -  mass  ■  is. 

•  _A(ii-\  '  -  mass  o'  i  r-  .  '  •  -  water 
.  sAx  -  mas.-;  of  ci.-'-pl.n  ...  sediment 
aecelernt  ion  of  gr  .  .’  >■ 

Assuming  all  coefficients  •  •.w-t.-inc,  the  equation  of  motion  is 
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wiie re  >:j'  •-  M  -l‘;v  :  tee  ,..e.  •'  t mi's  and  the  sediment  and  water  that 

ir.ovi  wi  ti:  it  in.vo.is  in»;  :t  inertia.  A  to  net  know  how  to  evaluate  My  at 
t  a  is  t  ir.a  t  see  ?x  :  I'scv  it.-  ..  *.  i  i'rot  i  .s  • )  >  and  '.ill  carry  it  through 

the  der  i  v;>  t  i  e:  . 
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>.  =  (>i,-e  Ah)g  /  iK-t  (i  -  )A  i  .  (4) 

t  i  w  s  w  g 

Let  \  =■  ,  i  |s  •  deflect  i  on  of  tile  OP.S  from  equilibrium  with  respect 

to  the  ocean  hot  to-m.  !>■  variable  v  is  thus  the  error  signal  between  true 
motion  oi  the  ocean  hett'.r.;  and  motion  of  tile  OBS . 

Substituting  v  tor  \  in  Kq.  iD, 
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The  term  on  the  right  containing  y  will  certainly  be  insignificant 
compared  to  the  others;  thus. 
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where  M  and  M  are  the  average  displaced  water  and  sediment  masses, 
w  s 

hq.  (6)  can  be  rewritten: 

? 

v+2h  u  v-Hi'  “v  =  Cz  (7) 

c  C '  C  '  s 

where  h  =  1/2Q  ,  the  damping  coefficient  of  the  bottom  coupling. 
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is  the  natural  angular 
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frequency  of  the  bottom  coupling)  and 


M,-M  -M 
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C  =  - -  ,  the  coupling  coefficient. 


Taking  the  transform  of  (7) 
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The  position  of  the  bottom  of  the  OBS, 
z ,  =  z  -x  =  z  - v - x 

Therefore  the  input  sicnal  to  tiie  OBS  (motion  of  the  case)  is 


i  =  z  t  +x 


z  -v 
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whose  transform  is 


i  =  ?.  -Y  , 


and  tin  transfer  functions  of  the  OBS  package  is 
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From  Eq.  (11)  it  is  evident  that,  given  the  constants  C,  hc,  and  uic , 
the  response  is  completely  described.  The  amplitude  response  curves  for 
six  cases  of  Eq.  (11),  shown  in  Figures  E-8  and  E-9,  demonstrate  good 
response  when  ui<uic ;  hc"i;  and/or  CO.  The  coupling  coefficient,  C,  can 
he  determined  with  the  aid  of  Eq.  (9)  provided  a  reliable  value  for  M* 
can  be  obtained  (see  Zelikovita  and  Prothero,  1980).  The  coupling 
frequency,  u>c ,  can  be  estimated  by  using  Eq .  (8).  Experimental  values  for 
. c  and  hc  can  be  obtained  from  mechanical  transient  tests  to  be  described 
in  the  following  section. 

When  conditions  are  such  that  C=0,  coupling  is  perfect,  and  the  OBS 
moves  exactlv  with  the  ocean  bottom,  although  the  package  may  still  have 
a  resonance  ai  uic  for  other  types  of  input  (such  as  the  transient  tests) . 
0=0  is  the  neutral  buoyancy  case  where  the  mass  of  the  displaced  water 
and  sediment  equal  the  mass  of  the  OBS.  Earlier  workers  studying  coupling 
in  swamps  ignored  the  effects  of  buoyancy  and  the  sediments  that  move  with 
the  instrument,  thus  arriving  at  the  case  where  C=l.  We  see  from  Eq.  (9) 
that,  since  M*  will  always  be  greater  than  Mj  and  since  Mw  and  Ms  will 
probably  never  be  negligible,  C  will  never  reach  a  value  of  1.  Note  that 
Eqs.  (8)  and  (9)  do  not  require  that  the  OBS  be  cylindrical  in  shape  to 
be  valid.  As  long  as  the  area  in  contact  with  the  bottom  does  not  change 
with  motion  of  the  OBS,  both  equations  will  hold. 

The  buoyancy  effect,  in  addition  to  modifying  the  shape  of  the 
response,  modifies  the  resonant  frequency  and  damping  from  that  obtained 
for  the  "swamp”  case  by  previous  authors.  From  Eq.  (8)  we  see  that  the 
resonant  frequency  is  the  RMS  of  that  produced  by  the  effective  spring 
constant  of  the  bottom  alone,  and  by  the  buoyancy  effect  alone.  As  a 
numerical  example,  for  the  HIC.R  geophone  package  (Sutton  et  al.,  1980)  in 
average  sediments  the  increase  in  nc  would  appear  to  be  less  than  1%. 
However,  the  effect  could  be  important  for  systems  with  large  surface  areas 
in  contact  with  the  bottom. 


We  can  investigate,  further,  the  effects  of  bearing  area  on  the 
coupling  frequency  by  rewriting  Fq.  (8). 
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r  represents  bearing  area 

r ' 1  represents  the  volume  <>:  sediment  and  water  entrained  with 
the  OBS 

c  [  depends  upon  tin  n'nl  i  gum  t  ion  of  the  OBS 
Mj  is  OBS  mass 

as  is  proportional  to  sediment  rigidity  modulus 
bs  depend,  upon  sediment  density  . 


where 
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E-8.  OBS  package  seismic  response  from  Eq  .  (11)  for  underdamped  coupling 
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PACKAGE  RESPONSE 

F.-9.  OBS  package,  seismic  response  from  Eq .  (11)  for  critically  damped 
oupl i ng . 
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2  3 

The  r  terra  results  from  buoyancy  and  the  r'  term  from  the 

entrained  sediment  and  water.  Assuming  r*  is  proportional  to  r,  for  a 
given  set  of  constants  (including  instrument  mass),  a>c  varies  as  rl/2 
and  r_l/2  for  small  and  large  r,  respectively,  and  varies  at  most  by 
either  r  or  r~3,  depending  upon  the  values  of  the  constants,  for  inter¬ 
mediate  values  of  r. 


file  case  C=1  mentioned  earlier  is  appropriate  for  the  horizontal 
response  of  the  OBS  package  to  horizontal  input  signals  (provided  by  the 
rigidity  of  the  bottom,  onlv).  For  horizontal  input.  Eq .  (11)  reduces  to 
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(12) 


In  this  case  the  response  decreases  continuously  at  frequencies  higher 
than  the  coupling  frequency;  as  !/*.•  for  frequencies  above  v)I/2hc. 


Mechanical  Transient  Tests 


A  schematic  diagram  for  a  simple  mechanical  transient  test  to 
determine  the  frequency  and  damping  of  OBS-bottom  coupling  (u>c  and  hc)  is 
•shown  in  figure  H-10.  This  test  is  analogous  to  the  classical  weight 
lift  calibration  for  earthquake  so t smographs .  The  float,  which  provides 
a  constant  upward  or  sideward  force  on  the  OBS  package,  is  released  at 
t.=0  when  the  current  through  an  electromagnet  is  interrupted.  With  use 
of  a  small  magnet  armature  and  light  string  or  thread  the  step  in  force 
can  be  made  quite  pure.  However,  the  pulley  rod  used  for  the  horizontal 
test  can  generate  a  second  signal  when  the  armature  hits  the  pulley.  In 
soft  sediment  a  separation  of  the  pulley  from  the  OBS  of  1  or  (preferably) 
2  m  is  generally  adequate  to  separate  the  primary  and  secondary  signals. 
The  size  of  the  float  can  be  adjusted  to  provide  a  signal  in  the  normal 
operating  range  of  the  OBS  and  to  test  for  possible  system  non- 
1  inear i t ies . 


An  upward  or  sideward  step  of  force  on  an  OBS  at  t=0  produces  initial 
conditions  of  z^(0)  =  zjq,  z^(0)  =  0  and  zj(0+)  =  Kzjq/M  for  the  position 
of  OBS  base  relative  to  its  equilibrium  position  at  rest  on  the  sediment. 
We  assume  that  the  OBS  moves  as  a  rigid  mass,  that  the  motion  is  recti¬ 
linear,  and  that  the  motion  of  tiie  inertial  masses  in  the  geophones  does 
not  affect  the  motion  of  the  OBS  (i.e.,  no  back  coupling). 

The  equation  of  motion  for  the  OBS  is: 


M,v,,viVi 


0  (assuming  no  seismic  input) 


(13) 


Mf  =  OBS  mass  plus  mass  of  sediment  and  water  included  in  the 
motion  of  the  OBS 


where 


'it;.  E-10.  Schematic  for  mechanical  step  of  force  test  to  determine 


“'c  and  hc' 
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D  =  damping  constant 

Kg  =  effective  spring  constant  of  the  OBS-bottom  coupling, 
including  any  buoyancy  effects. 

Rewriting  (13)  we  have 
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where  oc  =  Ks/Mj  and  hc  =  1  / 2 Q c  =  damping  coefficient,^  is  the  free 
angular  frequency  of  the  bottom  coupling,  and  Qc  =  l/2hc  gives  the  ampli¬ 
fication  (and  bandwidth)  of  any  coupling  resonance. 


Solutions  for  our  initial  conditions  are: 
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The  equation  of  motion  for  the  seismometer  mass  is: 
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where  zjq  is  the  relative  displacement  between  inertial  mass  and  OBS  frame. 
Taking  the  Laplace  transform  of  (19)  and  (for  example)  assuming  critical 
damping  hj.j  =  1,  the  seismometer  transfer  function  is: 
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The  emf  generated  in  the  seismometer  coil  is 
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and  the  signal  transfer  function  is 
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From  Eqs.  (18)  and  (18)  the  transforms 
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of  the  transient  input  signal  are 
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Combining  (23)  or  (2k)  with  (22)  we  obtain 
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,  cii>nals  for  the  mechanical  transients 

We  obtain  the  geophone  output  sign 

by  taking  the  inverse  transform  of  (2o)  and  (26). 
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Figure  E- 1 1  shows  theoretical  signals  from  mechanical  transients  for 
live  different  cases  listed  in  Table  E-l, and  Figure  E-12  shows  amplitude 
spettra  for  three  of  the  cases.  These  results  are  appropriate  for  hori¬ 
zontal  .and  vertical  inputs  recorded  on  horizontal  and  vertical  geophones, 
respectively.  In  addition,  bv  observing  the  vertical  geophone  output  from 
horizontal  transients  and  vice  versa  the  susceptibility  of  an  OBS  to  cross¬ 
coupling  distortion  can  be  evaluated.  Some  results  from  mechanical 
transient  tests  conducted  in  shallow  water  off  Oahu  were  presented  by  Sutton 
et  aL.,  (1977b  and  1978).  Sample  results  from  the  Lopez  Island  experiment 
f o 1 1 ow . 


lahle  E-l.  OBS  coupling  parameters  for 
curves  shown  in  Figures  E-l 1  and  E-12. 
alui  hM  'lre  l'>oth  ^Oinl  to  1.0. 
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A  great  deal  ot  information  on  coupling  is  contained  in  the  digital 
records  from  the  mechanical  transients  and  from  the  blasting  cap  and  airgun 
shots  obtained  during  the  l.opez  experiment.  These  results  are  summarized 
in  Sutton  et  a!.,  (1980)  and  discussed  further  in  Zelikovitz  and  Prothero 
(  1980),  I  ii  I  hi  1  1  et  al.  (1980),  Lewis  and  Tuthill  (1980),  and  .Johnson  and 
McAlister  ( 10801. 

In  general  ,  l  he  responses  te  t  Ik  transient  tests  ot  the*  instruments 
special  iv  designed  lor  use  as  standards  at  l.opi  Island  exhibit  broad, 
rather  featureless  spectra,  whereas  most  of  the  spectra  from  tests  on  the 
operational  (IBS's  show  one  or  more  spectral  peaks  with  various  bandwidths 
and  .amp]  ifiration  (Fig.  1.-13).  Most  of  the  OBS's  show  a  peak  between  about 
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Fig.  K—  II-  I'lii’ii re t  i ca  1  gcophone  outputs  for  mechanical  transient  test  cases 
one  through  seven  listed  in  Table  F-l  . 
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.  F.-I2.  Spectra  for  three  of  the  test  outputs  (2,4,  and  6)  shown  in 
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three  are  >  perational  OBS's 
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X  i n cl  1J  Hz  tor  both  the  vertical  and  horizontal  inputs,  suggesting  some 
kind  ol  geological  wave- guide  resonance;  however,  its  absence  on  the 
special  standard  instruments  and  genera]  agreement  with  the  coupling  theory 
argue  that  it  is  a  coupling  resonance  and  that  the  DBS's  have  roughly 
equivalent,  hot  tom-coup  1 i ng  design  parameters.  In  this  respect  the  MITEV 
and  HIDBY  (both  having  geophones  separated  from  the  main  package)  appear 
to  he  closest  in  response  to  the  standard  instruments. 

Comparison  of  the1  spectra  from  the  mechanical  transient  tests  with 
those1  Iron  cap  and  airgun  shots  and  from  the  background  noise  often  show 
similar  peaks  and  attenuation  oi  higher  frequencies,  in  agreement  with  the 
coupling  theorv  (Figs.  F.-H  and  t-9) .  The  Q  ol  the  coupling-resonances  are 
generally  greater  than  predicted  from  the  published  theories  (Fig.  F-l) 
that  do  not  consider  buoyancy  forces.  The  coup  1 ing- resonance  frequencies 
mnv  he  more  reliable  predictable. 

Figuri  F.-l'-t  shows  the  signals  recorded  digitally  from  an  airgun  shot, 
approximate! v  HO  m  southwest  of  the  center  of  the  DBS  array.  Considering 
differences  in  sensitivity  and  frequency  response  of  the  various  OBS’s  the 
first  few  cedes  of  the  signal  recorded  on  the  vertical  components  agree 
well,  although  some  arrivals  seem  somewhat  more  emergent  than  others. 
Between  the  first  arrivals  and  the  large  "granddaddy"  waves  on  the  right,  a 
wave-train  with  predominantly  horizontal  (mainly  north-south)  motion  is 
observed  on  the  plate  (PH  and  neutral  density  (NO)  standards  and  on  HIGBH. 
Although,  the  signal  is  hardly  observable  on  the  vertical  standards,  it  is 
clearly  recorded  on  nest  a;  the  DBS  vwrt  icals.  This  signal  appears  to  be 
an  example  of  cross  coupling  between  a  horizontal  input  and  vertical  output. 
The  observed  signal  ■  to*-  i  g  i  ■  .  u  in  t  i  intent  are  quite  dependent  upon  the 
stiffness  of  the  Shea  I  ike.  sediments  and  upon  the  frequency  of  tile  input 
signal  (narrow  hand  near  i*  Hz'.  For  example,  BIDV  exhibited  a  strong 
response  near  *>  11/,  to  a  iior  i  zontal  transient  input  (Fig.  1:1— 15),  close  to 
the  principal  frequency  of  the  airgun  signal.  Other  instruments  that  show 
less  relative  cross  coupling  I  mm  the  airgun,  e.g.,  I'CSBV,  also  shown  in 
Figure  E-lb,  exhibited  equal lv  strong  responses  on  the  vertical  sensors  to 
horizontal  transients,  but  at  different  frequencies.  Lewis  and  Tuthill 
(H801  present  a  more  complete  discussion  of  the  cross-coupling  phenomenon. 


(lone  !  1 1 s  i  on 


The  theorv  and  mechanical  transient  tests  outlined  here  and  the  results 
of  the  Lope::  Island  f  nt  <•reomp.tr  i  sen  experiment  provide  reliable  guidelines 
tor  p r id i v  t  ing  tlu  n sponsi  of  an  DBS  under  various  bottom  conditions. 

This  kn--wl edge  makes  it  possible  to  establish  parameters  for  ocean  bottom 
seismometer  design.  < toner. i! !y,  a  low  smooth  prof  ile  in  the  water,  a  large 
bearing  are  t,  and  a  low  package  density  should  provide  optimum  coupl  ing 
to  tin  ha*  t  >:•:  in  soft  sediment  .  Careful  design  should  result  in  a  package 
that  will  couple  weM  to  anv  t  v|k  of  ocean  bottom. 
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H- >>•  i  znni  .1 1  mocltanic.il  transient  tests  recorded  on  vertical 
s.  iloriz.m'a]  line  under  t  i  ni<e  signals  represents  two  seconds 
I  lines  l> racket  portion  mily/.ed.  Spectra  are  in  dB  and  Hz; 
■spectra  are  for  a  noise  sample  of  tlie  same  length  as  signal, 
i  i nned  i . 1 1 . ■  1  v  pro  eding  transit’ll l  . 
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Introduc  t ion 


In  seismology,  it  is  desired  to  record  ground  motion  accurately. 
However,  Ocean  Bottom  Seismometers  (CBS)  ,  which  are  usually  rather  massive 
instruments,  can  significantly  distort  the  measured  ground  motion, 
particularly  when  they  rest  on  sediments  of  relatively  weak  stiffnesses. 
This  response  problem,  which  has  been  recognized  in  the  case  of  land 
seismometers  for  over  30  years,  is  known  to  depend  on  the  physical  pro¬ 
perties  of  the  instrument  and  of  the  soil  on  which  the  instrument  rests. 

It  can  be  more  acute  for  an  OBS  due  to  the  possibility  of  deployment  on 
softer  underlying  sediments,  and  more  complex  due  to  water-instrument 
interactions . 


The  original  theoretical  development  in  soil-structure  interaction 
(Kichart  et  al. ,  1970)  provides  a  framework  within  which  we  can  study  OBS 
response.  Work  by  K.  L.  Hamilton  (1970,  1971a, 1971b)  suggests  treatment 
of  marine  sediments  as  an  elastic  solid  that  can  be  represented  by  a 
spring-mass-damper  svstem.  (Normally,  soil-structure  interaction  studies 
require  the  use  of  finite  element  techniques  because  they  most  often 
involve  structures  with  many  degrees  of  freedom, resting  on  complex 
foundations.)  The  stiffness  and  damping  coefficients,  or  foundation 
impedances,  arc  functions  of  the  bearing  radius  of  the  instrument  and  the 
elastic  properties  of  the  soil.  The  presence  of  water  can  be  accounted 
for  by  using  an  "effective  mass”  for  the  OBS,  which  includes  the  inertial 
effects  of  the  surrounding  fluid.  The  theory  and  data  analysis  deal 
strictly  with  the  vertical  motion.  Horizontal  motions  will  also  produce 
important  effects  hut  art'  not  treated  in  tiiis  paper. 

Since  the  ultimate  goal  of  this  study  is  finding  a  calibration 
technique  for  any  OBS  at  any  site,  we  explore  the  response  problem  both 
through  modeling  experiments  and  by  analvsis  of  the  Lopez  Island  vertical 
transient  tests.  The  Lopez  results  allow  us  to  compare  the  output  of 
many  instruments  and  standards  subjected  to  identical  step-function  inputs. 

The  modeling  experiments  similarly  consider  the  response  of  an  ideal 
model  to  a  variable  frequency  sinusoidal  input.  Foam  rubber  of  known 
elastic  properties  serves  as  a  homogeneous  "ideal"  soil  in  the  modeling 
experiments.  We  investigate  so i 1— instrument  interaction  bv  applying 
vertically  oscillating  forces  to  masses  of  varying  bearing  radius  placed 
on  the  foam  soil  model.  This  enables  us  to  verily  the  continuum 
mechanical  theory  lor  the  simple  eases  and  model  more  complex  configura¬ 
tions,  which  do  not  easily  lend  themselves  to  a  theoretical  treatment. 

The  experiments  with  foam  rubber  also  help  us  to  identify  some  of  the 
important  considerations  in  the  analysis  ol  the  Lopez  data.  In  addition, 
they  suggest  that  a  possible  in-si tu  technique  for  calibrating  an  OBS 
might  make  use  of  an  attached  shaking  unit. 

A  harmonic  oscillator  subjected  to  a  step-function  input  force  will 
exhibit  a  decaying  sinusoidal  output.  Clearly,  the  Lopez  transient  test 
data  shows  this  decaving  oscillatory  behavior  with  a  measurable  frequency 
tor  every  instrument.  The  decay  rate  can  also  be  measured  and  is 
function  of  the  damping  in  the  system.  By  relating  these  ringdown 
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frequencies  to  the  effective  masses  and  bearing  radii  of  the  OBS's,  a 
shear  modulus  value  for  Shoal  Bay  mud  can  be  obtained.  Since  the  OBS's 
and  standards  were  placed  on  the  bottom  in  a  circle  of  small  diameter, 
the  shear  modulus  value  felt  by  each  instrument  should  be  comparable. 

We  attempt  to  explain  four  notable  discrepancies  by  examining  the  assump¬ 
tions  inherent  in  the  simple  model.  Embedment  due  to  bearing  pressure 
of  the  instrument  in  a  nonuniform  vertical  soil  profile  provides  a 
convincing  explanation  of  why  certain  instruments  feel  a  higher  shear 
modulus . 

Damping  factors  can  also  be  obtained  from  the  Lopez  data.  Although 
the  logarithmic  decrement  of  the  ringdown  cannot  be  measured  with  high 
accuracy,  a  range  in  value  for  each  damping  factor  is  reported.  These 
values  are  compared  to  a  theoretical  damping  curve  for  Lysmer's  and 
Kichart's  (196b)  elastic  half  space  analog. 


Background  and  Theory 

K.  i..  Hamilton  in  several  papers  (1970,  197.1a,  1971b)  concludes 
that  equations  ol  Hookean  elasticity  can  be  used  to  calculate  the  elastic 
constants  for  water-saturated  sediments.  Therefore,  we  begin  by  model¬ 
ing  an  i IBS  as  a  rigid  circular  mass  resting  on  an  elastic  half  space  with 
properties,  L-shear  modulus,  v-Po  isson ’ s  ratio,  and  n-density  (see  Fig.  F-l). 
Hsieh  (1962)  developed  the  lumped  parameter  system  shown  in  Figure  F-l  as 
f o I  Lows : 

Mz  -()-!’  (1) 

i 1  =  C  +  K  X  ( 2 ) 

z  z 

Mz  +  C  7.  +  K  7i  =  0  =  q„e1Wt  (3) 

z  z  '  0 


Lysmer  and  Richart  (1966)  found  frequency  independent  vertical  stiffness 
and  damping  constants  as  functions  of  the  mass  radius  (R)  and  soil 
properties  (i,  ,  and  ,>.  The  constants  Kz  and  Cz  are  given  by: 


4(1  K  =  3 . 4R 

1-v  ’  z  T-V 


(4) 


Substituting  these  constants  into  Eq -  (3)  yields  a  linear  differential 
equation  witli  constant  coefficients,  sometimes  called  Lysmer's  analog. 

Using  Eq. (4),  the  natural  frequency  of  this  damped  harmonic  oscillator  will  be 


(3) 

where  D  is  a  damping  constant  which  is  the  actual  damping  constant  divided 
by  the  critical  damping  constant  C’^./C^.  Lysmer's  analog  breaks  down 
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for  D  >  .425.  An  excellent  discussion  of  this  subject  is  found  in 
Richart  et  al.  (1970). 

To  this  point,  the  instrument  water  interaction  has  been  neglected. 
There  are  three  terms  that  enter  the  equation  of  motion:  hydrodynamic 
drag,  wavemaking,  and  virtual  mass.  The  size  and  depth  location  of  the 
instrument  make  the  relative  magnitude  of  drag  and  wavemaking  forces 
insignificant  with  respect  to  the  inertial  effect  of  the  fluid  (Byrd,  1978). 

Figure  F-2  shows  the  inertial  and  mechanical  forces  acting  on  a 
submerged  mass  resting  on  an  elastic  solid,  which  is  represented  by  a 
spring  and  dashpot.  The  position  of  the  mass  relative  to  an  inertial 
frame  of  reference  is  x  +  Rq  +  AR,  where  Rq  is  the  steady  state  radius 
of  the  earth,  AR  is  the  free  field  displacement,  and  x  is  the  displace¬ 
ment  of  the  mass  relative  to  the  free  field  position  of  the  seafloor. 

The  acceleration  of  M  relative  to  the  seafloor  is  then  x.  There  is  also 
a  virtual  mass  term  or  added  mass  effect  from  the  movement  of  the  mass 
through  water.  This  effect  is  caused  by  hydrodynamic  forces  exerted  on 
an  accelerating  mass  because  of  the  necessary  acceleration  of  the 
surrounding  water  through  which  the  mass  moves  (Batchelor,  1967).  Since 
the  wavelength  of  sound  at  100  Hz  is  about  15  m,  we  can  assume,  for  OBS- 
sized  instruments,  that  the  water  surrounding  the  instrument  moves  with 
the  bottom.  So,  the  virtual  mass  force  due  to  the  relative  motion  of  the 
instrument  and  water  will  be  Mvx.  Now  the  buoyancy  force  is  due  to  the 
difference  in  the  weight  of  the  OBS  and  the  weight  of  the  displaced  water, 
which  is  its  mass  multiplied  by  the  acceleration  of  its  local  frame  of 
reference.  We  must  include  the  acceleration  of  gravity  (g)  and  the 
acceleration  because  of  ground  motion  AR ,  so  the  total  buoyant  force 
is  given  by  (M  -  Mw)(g  +  AR). 

The  use  of  this  buoyancy  term  is  the  simplest  way  to  include  inertial 
effects  of  the  acceleration  of  the  water-OBS  system.  The  equation  of 
motion  obtained  by  balancing  these  forces  and  the  stiffness  and  damping 
forces  is: 

Mx  +  Mx  +  Cx-K  (y  -  x)  +  (M  -  M  )  (g  +  AR)  =  0  (6) 

V  Z  7.  U  W 

where  M  is  the  mass  of  the  instrument,  Mv  is  the  virtual  mass  effect  due 
to  the  acceleration  of  surrounding  water.  My  is  the  mass  of  the  displaced 
water,  xo  is  the  equilibrium  position  of  the  spring,  and  A’k  is  the  free 
field  ground  motion  acceleration.  The  weight  of  the  instrument  in  water  is 
balanced  by  the  equilibrium  spring  force.  Therefore  the  equation  reduces  to 

(M+M)x+Cx+Kx=  -Ah(M  -  M  )  .  (7) 

V  z  z  w 

We  see  that  if  the  instrument  is  neutrally  buoyant  (M  -  =  0) , 

x  and  x  will  be  zero  as  there  is  no  input  force  to  the  mass  spring  system. 
Thus,  the  OBS  moves  exactly  with  the  bottom  as  expected  intuitively. 

A  similar  equation  of  motion  has  been  derived  by  Sutton  et  al.  (1980, 
Eq.  (5)).  Our  equations  differ  in  that  we  include  virtual  mass  effects 
from  hvdrodvnamic  theory  and  ignore  their  term  Ms,  which  is  the  mass  of 


tin'  displaced  sediment,  which  i>ivos  rise  to  an  additional  Xuovancy  foive. 
1'his  force  may  ho  neglected  for  the  Lopez  experiment.  These  oq  tat  inns  of 
motion  lead  to  a  roupl  inn  constant,  C  =  (M  -  M._()  (M  Mvi,  whim 
represents  the  ratio  of  the  buoyant  mass  to  the  dynamic  mass  of  the  OBS- 
water  system  ("Sutton  et  al.  19R0a).  C  ~  0  represents  nerfect  coupling, 
and  (1  =  1  is  the  worst  ease.  Stated  another  way,  the  op.s  rosponst  ran  be 
compared  to  that  of  a  harmonic  oscillator  _th  input  force  proportional 
Co  the  unbalanced  huovancv  of  the  OHS  in  the  water  ami  a  mass  equal  to 
tile  DBS-water  dynamic  mass.  The  spring  constant  and  damping  are  caused 
by  the  soil  stiffness  arid  energy  radiated  to  infinity. 

Since  the  free  field  ground  motion  (.'  K)  is  zero  for  the  vortical 
transient  teats  at  Lopez,  the  equation  is  modified  by  replacing  the 
right-hand  side  by  the  force  induced  by  a  step  function  transient  input: 

( M  +  M  ) x  t  C  X  +KX~l(t).  ( 3 ) 

V  7.  Z 

Tor  sinusoidal  oscillation  of  the  instrument  with,  an  attached 
shaking  unit: 

(M  M  )  x  +  C  X  +  K  X  =  -M  v  (9) 

t'  V  Z  Z  s' 

whore  M.;  -  mass  of  shaking  unit 

v  is  the  aocelerat ion  of  the  shaking  mass. 

!>;.  (/>)  gives  the  natural  f recpicncv  for  anv  of  the  above  three  cases  if  M 
is  replaced  bv  the  e! foctive  mass  Mr  where 

M_  =  M  ¥  M  .  (10) 

R  v 

Those  I amil iar  with  various  OBS's  realize  that  they  are  generally 
complex  in  form;  uowvvcr,  most  can  he  approximated  by  combinations  ot 
simple  geometric,  shapes.  The  added  mass  terms  for  a  sphere  O/2  M,0)  »  a 
disc  (4  ,  ,WK * ) ,  and  n  cvl-'nder  (,6MW)  are  used  to  calculate  the  effective 
mass  oi  each  instrument  and  standard  (Batchelor,  1967 ;  Byrd,  1973).  Mw 
is  the  mass  oi.  the  displaced  water,  .  w  is  the  water  density,  and  R  is  the 
rail  i .is  of  t !  e  :1  iso. 

L.-s.  (?)  ,  (3),  and  (9)  assume  a  mass  with  a  circular  hearing  radius. 

Tor  vertical  motion,  the  theory  is  valid  for  all  r ripod-oased  instruments 
»,•  1 1  h  1  intiiar  contact-.  As  well,  OBS's  where  one  dimension  does  not  exceed 
twice-  ;  iii-  other  r.iv  be  similarly  treated  (Ri  chart  ct  al.,  1970).  (The 
bearing  radius  is  computed  bv  assuming  the  rectangular  area  to  he  circular.) 
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Damping  in  Lvsmer's  analog  is  due  to  radiation  losses  by  the 
propagation  of  seismic  energy  to  infinity  and  can  be  calculated: 


C 


c 


(ID 


2VK  M  =  2/4GRM 
z  h  /  1. 

/  1-v" 


(12) 


!)  =  .42  5  /  v  11 


(13) 


where  Cc  is  critical  damping  and  B  is  a  dimensionless  mass  factor 


B 


1-. 


4 


(14) 


B  is  a  comparison  of  the  mass  of  the  instrument  to  the  mass  of  the 
soil  influencing  the  instrument.  In  terms  of  a  harmonic  oscillator  model 
it  is  a  ratio  of  the  moving  mass  to  the  mass  of  the  spring.  A  high  B 
value  signifies  a  comparatively  massless  spring  and  therefore,  correspond¬ 
ingly  low  radiation  damping.  However,  low  B  values  signify  a  spring  with 
mass  approaching  that  of  the  instrument,  and  thus  high  radiation  damping 
would  result.  The  analogy  to  a  mass  sprang  system  differs,  however, 
since  the  energy  is  radiated  away  to  infinity  in  the  half-space  but  it  is 
contained  in  modes  of  oscillation  of  the  spring  in  the  mass-spring  system 
The  effect  of  the  mass  of  the  half-space  that  moves  with  the  OBS  is 
included  in  F.q.  (3)  in  the  damping  coefficient,  i).  It  should  be  noted 
that  under  the  conditions  that  hysmcr's  analog  holds  (B  >  1),  the  mass  of 
the  i IB S  must  he  relatively  large  compared  to  Lite  mass  of  the  sediment 
that  would  move  under  its  in’.iuence. 


The  above  theory  includes  several  idealizations.  OBS's  do  not  rest 
on  the  bottom  sediments  hut  embed  themselves  to  a  depth  depending  most 
importantly  on  their  hearing,  pressure.  Also,  the  soil  beneath  the 
instrument  is  neither  homogeneous,  isotropic,  nor  semi-inl xnite .  These 
assumptions  and  others;  of  a  more  subtle  nature  will  be  discussed  further 
in  tiie  analysis  and  discussion  section. 


Model_Invest  ignti on 

For  i  ■.  ,.mer '  s  and  Kichart's  (196b)  mass-spring-damper  model,  a 
reciprocity  exists  between  response  due  to  ground  motion  (from  l.q.  7): 


M-M 


M+M 

v 


"0 


! 


1  +  i 


z 

(M+M  , )  . 

V 


(15) 
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and  the  response  computed  bv  shaking  the  instrument  with  a  displacement 
y  =  v  exp (i. it): 


The  only  dillerence  is  an  amplitude  factor  reflecting  the  inertia  transfer 
between  the  shaker  and  the  larger  instrument  mass. 


In  order  to  verify  the  theorv  and  develop  a  procedure  for  studying 
more  complex  cases,  experiments  were  performed  with  foam  rubber.  The 
experiments  were  conducted  in  air  (My,  M^.  =  0)  and  were  intended  to  study 
soil-structure  interaction.  Different  masses  of  varying  bearing  radius 
were  placed  on  foam  rubber  and  oscillated  with  a  sinusoidal  input  force. 
The  resonant  frequency  depends  on  the  stiffness  and  damping  constants  and 
therefore  on  the  elastic  properties  of  the  foam  (Eq.  4).  Young's  Modulus, 
the  modulus  of  rigidity  and  the  density  of  several  types  of  foam  were  also 
measured  in  static  experiments.  A  high-density  loaded  foam 
(Cl  =  1.5  x  10’N/M‘  ,  E  =  1.1  x  10  N'/M'  ,  and  .  =  11  kg/M')  was  chosen  for 
use  in  all  experiments  because  of  its  linear  stress-strain  relationship 
over  a  wide  range  of  stresses.  We  then  compared  values  of  the  elastic 
parameters  determined  from  the  static  tests  to  those  determined  by  shak¬ 
ing  a  model  OHS  resting  on  foam. 


Kigurt  F-l  shows  the  experimental  configuration  used  in  finding  the 
resonant  t  requeues  of  the  mass-spring-damper  analog.  A  sinusoidal 
voltage  is  applied  to  geophone  A,  causing  the  moving  mass  to  oscillate 
and  tnus  applying  a  sinusoidal  force  to  the  system.  f.eophone  15  acts  as 
the  seismometer,  which  senses  the  motion  induced  by  Geophone  A.  Both 
geophones  are  mounted  rigidly  in  a  PVC  holder  that  is  fastened  to  discs 
of  various  radii.  The  mass  of  the  model  includes  the  mass  of  the  geo¬ 
phones,  holder,  disc,  and  any  added  weights.  The  area  of  the  foam,  on 
which  the  model  rests,  is  large  enough  to  approximate  the  elastic  haif 
space  in  two  dimensions.  The  thickness  of  the  foam  can  be  varied  to 
allow  us  to  model  the  transition  between  an  elastic  half  space  and  an 
elastic  layer  on  a  rigid  stratum.  By  slowly  varying  the  input  frequency, 
the  resonant  frequency  is  determined  for  each  disc-mass  system.  Plotting 
the  natural  frequence  against  (R/M)’/2  should  return  a  lino  proportional  to 
the  dynamic  shear  modulus  of  the  foam  according  to  Eq .  (5),  assuming 
small  damping.  The  results  are  shown  in  Figure  F-4  along  with  a  hashed 
line  corresponding  to  the  shear  modulus  obtained  in  static  tests.  One 
source  oi  variation  in  the  points  from  the  static  shear  modulus  value  is 
tin  neglecting  damping.  The  radiation  or  geometric  damping  computed 
using  Eqs.  (b)  and  (10)  is  negl  gible  because  of  the  low  density  of  the 
foam;  however,  the  viscous  and  material  damping,  although  unaccounted  lor 
in  ,in  elastic  model,  will  cans-  deviations  of  the  points  from  the  line. 

Tiie  ratio  of  the  thickness  of  the  foam  (H)  to  the  bearing  radius  of 
the  disc  (K)  is  also  tabulated  for  each  point  on  Figure  F-4.  For  H/R 
large,  the  loam  can  he  treated  as  a  semi- inf  ini te  elastic  half  space,  but 
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F-3.  Schematic  of  a  foam  rubber  modeling  experiment,  designed  to  investigate  instrument  soil 
interactions.  Geophone  A  is  driven  by  an  external  sinew nve  current  source  and  Geophone  B  measures 
the  resulting  displacement. 
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for  large  bearing  radius  discs,  the  rigid  stratum  beneath  the  foam  must 
be  taken  into  account.  Note  points  2,  3,  and  6  which  show  progressively 
larger  shear  modulus  values  with  respectively  decreasing  H/R  values. 


Lopez  Data  Interpretation 

The  elastic  soil  model  was  also  applied  to  the  Lopez  vertical 
transient  test  data.  The  bearing  radius  and  mass  of  each  instrument  was 
calculated  from  the  physical  properties  supplied  by  its  institute  of 
origin  (Table  F-l  ).  The  error  estimate  associated  with  each  OBS  effec¬ 
tive  mass  is  based  on  the  degree  to  which  it  could  be  thought  of  in  terms 
of  simple  geometric  shapes. 

The  natural  frequency  of  each  instrument-soil-water  system  was 
measured  directly  from  the  vertical  transient  test  data  acquired  at  Lopez. 
(For  data  playouts,  see  Loncarevic,  1979;  Sutton  et  al.,  1980b.)  Measured 
ringdown  frequencies,  as  recorded  in  Table  F-2, have  an  estimated  error  of 
.5  Hz  for  f £  |0  Hz  and  +  1  Hz  for  ffj  £  10  Hz.  f[j  is  plotted  against 
[R(l  -  2D2 ) /Mg]  in  Figure  F-5.  Two  lines  are  superimposed  to  illustrate 
the  trends  in  the  data,  including  approximate  shear  modulus  and  velocity 
values . 

Damping  factors  can  also  be  retrieved  from  the  original  transient 
test  data  by  modeling  the  decay  of  a  damped  harmonic  oscillator  (see  also 
Sutton  et  al.,  1980a).  The  decay  time  is  quite  sensitive  to  the  damping 
for  most  of  the  transients.  In  Figure  F-6  the  damping  factors  are  plotted 
against  B  values,  calculated  using  Eq .  (10)  and  the  data  in  Table  F-l.  For 
comparison  the  D  versus  B  curve  for  a  true  elastic  half  space  is  also 
plotted.  A  value  for  v  of  the  .48  is  used  throughout. 


Analysis  and  Discussion 

If  all  the  assumptions  implicit  in  the  proposed  model  were  valid, 
all  the  plotted  points  in  Figure  F-5  would  fall  on  a  single  straight  line 
corresponding  to  the  shear  modulus  of  Shoal  Bay  mud.  Also,  some 
variation  in  the  soil  beneath  each  instrument  is  expected  and  will  cause 
scatter  in  the  computed  shear  modulus  values. 

Up  to  now,  we  i  ave  treated  the  Shoal  Bay  mud  as  a  semi-infinite 
elastic  half  space.  We  do,  however,  have  reflection  data  (Sutton  et  al., 
1980b)  from  this  site  indicating  that  the  top  silt  layer  under  the 
instruments  is  only  21Ato  3  m  deep.  The  properties  of  the  "basement 
rock"  are  not  known.  Intuitively,  it  might  be  either  a  glacial  deposit  or 
a  hard  siltstone  or  sandstone,  but  it  is  definitely  of  higher  rigidity 
than  the  overlying  silt. 

Arnold  et  al.  (1955)  and  Bycroft  (1956)  studied  the  effects  of  a 
rigid  boundary  under  an  elastic  layer  and  Luco  (1974)  investigated 
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.  F-3.  Hiis  plot  of  ringdown  frequence  versus  [  K  <  1  -  JD  )/M£]  '  should  return 
the  shear  modulus  of  Shoal  Bay  mud.  The  two  lines  represent  the  shear  velocities 
of  13  and  34  m/sec,  respectively. 
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foundation  impedances  for  layered  elastic  media.  One  important  conclusion 
of  these  studies  is  that  the  underlying  layer  begins  to  have  a  signi- 
cant  effect  on  the  stiffness  when  H/R  is  less  than  6;  H  is  the  thickness 
of  the  top  elastic  layer  and  R  is  the  bearing  radius  of  the  instrument. 

This  effect  is  also  observed  in  the  foam  rubber  experiment. 

Another  important  effect  is  the  varying  embedment  of  the  instruments 
resulting  from  variations  in  their  bearing  pressure.  The  bearing  pressure 
in  Table  F-l  is  a  static  value  equal  to  the  weight  of  the  OBS  in  water 
divided  by  the  total  footpad  area.  If  the  instrument  were  deployed  in 
free  fall,  the  pressure  would  be  considerably  greater  when  it  hits  bottom. 
Considerable  variation  in  embedment  might  be  expected  due  to  this  effect. 
Figure  F-6  shows  a  plot  of  bearing  pressure  versus  shear  modulus,  G, 
inferred  from  the  ringdown  frequencies.  A  clear  trend  toward  increasing 
shear  modulus  for  more  massive  instruments  is  indicated. 

The  damping  factors  are  plotted  versus  B  in  Figure  F-7.  The  deviation 
from  the  elastic  half  space  model  increases  at  low  B  values,  where  Lysmer's 
analog  begins  to  break  down.  Sources  of  damping  such  as  hydrodynamic 
losses  and  losses  due  to  soil-water  exchange  should  be  considered.  It  is 
surprising  that  the  plate  standard  has  such  low  damping.  Its  very  low 
bearing  pressure  may  be  a  factor,  allowing  it  to  rest  on  the  very  weakest 
sediments,  yet  its  large  diameter  causes  it  to  feel  deeper,  more  compe¬ 
tent  material  for  the  transient  test.  Since  the  damping  is  due  to 
radiated  seismic  energy,  reflected  energy  from  underlying  layered  struc¬ 
tures  is  sure  to  cause  significant  effects.  Damping  factors  would 
intuitively  seem  to  be  more  sensitive  than  ringdown  frequency  to  reflected 
energy  from  layered  structures. 

With  geometric  damping  included,  the  average  shear  velocity  recorded 
by  the  instruments,  which  fall  along  line  "1"  in  Figure  F-5, is  15  m/sec; 
and  for  those  along  line  "2",  33  m/sec.  It  is  both  interesting  and 
encouraging  to  note  that  Lewis  and  Tuthill  (1980)  indicate  that  "granddaddy" 
waves  measured  at  Lopez  had  group  velocities  ranging  between  10  and  30  m/sec. 
Also,  the  time  difference  between  P-S  conversions  and  initial  P  arrivals 
from  airgun  records  suggest  a  shear  wave  velocity  of  15  m/sec  in  the  3-m 
silt  layer. 


In  Situ  Calibration 

The  reciprocity  between  shaking  the  instrument  and  ground  motion 
demonstrated  by  Eqs.  (15)  and  (16)  suggests  the  viability  of  an  in  situ 
calibration  technique.  A  first  step  might  be  to  verify  Eq.  (7)  by 
conducting  foam  rubber  modeling  experiments  in  water.  (These  experiments 
would  concentrate  on  the  instrument-water  interactions.) 

An  in  situ  calibration  method  might  rely  on  an  attached  shaking  unit 
to  determine  the  stiffness  and  damping  constants  of  the  soil,  following 
the  analysis  of  the  Lopez  data  and  foam  rubber  experiments.  An  estimate 
of  the  vertical  motion  response  function  would  be  computed  from  these 


elastic  half  space  analog 


Fig.  F-7.  The  damping  factors  measured  from  the  transient  test  data  are  plotted  against  the  effective 
mass  ratio  s  of  the  OBS's.  The  theoretical  curve  for  a  simple  elastic  model  is  illustrated  for 
comparison. 
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soil  constants,  the  virtual  mass  term  and  the  physical  properties  of  the 
instrument  and  Eq.  (15).  The  response  function  for  that  particular  site 
would  then  be  used  to  deconvolve  the  recorded  output.  More  work  is 
necessary  to  understand  the  instrument  response  to  horizontal  ground 
motion  inputs.  Instrument  rocking  and  shear  in  the  water  column  are 
expected  to  cause  significant  effects. 


Implications  for  OBS  Design 

Eq.  (7)  suggests  that  a  high  damping  factor  will  result  in  the  most 
smoothly  varying  and  uncomplicated  coupling  response.  Further,  if  Kz  is 
increased,  the  coupling  resonance  may  be  above  the  seismic  band  of  great¬ 
est  interest,  so  it  can  be  ignored.  This  would  imply  that  the  bearing 
area  should  be  large  to  decrease  the  B  value.  However,  the  plate  "feels" 
the  same  shear  modulus  as  do  the  other  instruments  yet  has  an  anomalously 
low  damping  factor.  Another  way  of  improving  the  coupling  is  to  make  the 
instrument  more  nearly  neutrally  buoyant.  This  will  affect  the  horizontal 
coupling,  although  probably  adversely.  Given  the  complications  in 
variable  soil  conditions,  impact  geometry  upon  deployment,  and  variable 
damping  factors,  an  in  situ  method  of  calibration  becomes  highly 
attractive  and  may  be  the  only  way  of  calibrating  an  OBS  tc  ground  motion. 


Conclus ions 


We  have  a  reasonable  understanding  of  the  Lopez  Island  vertical 
ringdown  frequencies.  Although  a  simple  uniform  elastic  model  does  not 
explain  all  the  results,  it  can  be  refined  to  give  qualitative  explana¬ 
tions  of  certain  discrepancies.  The  damping  factors  do  not  fit  the  half 
space  model  and  it  seems  that  reflected  radiation  or  other  factors  must 
also  be  important. 

The  model  fit  establishes  that  the  resonant  frequencies  excited  by 
the  Lopez  vertical  transient  tests  are  due  to  soil-instrument-water 
interactions.  These  resonances  occur  within  the  bandwidth  of  seismic 
interest  and  therefore  will  distort  the  recorded  signal. 

The  eventual  aim  in  studying  this  problem  is  to  develop  an  in  situ 
calibration  technique  for  any  OBS  at  any  site,  since  the  response  will 
surely  be  extremely  site  dependent.  The  transient  tests  performed  at 
Lopez  Island  and  the  foam  rubber  experiments  yii Id  the  response  of  an 
instrument  to  ground  motion.  All  of  the  dynamic  and  elastic  soil 
parameters  relevant  to  vertical  motion  are  contained  in  the  transient 
test  data.  In  situ  calibration  is,  therefore,  not  only  desirable  but 
feasible  and  should  be  considered  by  every  institute  deploying  OBS's. 

A  complete  understanding  of  OBS  response  also  requires  studies  of 
the  horizontal  motions  of  the  instruments.  This  is  a  more  complex 
problem,  yet  one  that  is  tractable.  Foam  rubber  modeling  would  be 
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particularly  useful  for  horizontal  motion,  to  demonstrate  the  behavior 
of  differing  instrument  bearing  configurations.  This  would  avoid  a 
lengthy  theoretical  development  that  is  not  critical  to  either  an  under¬ 
standing  of  the  problem  or  the  determination  of  the  appropriate  response 
deconvolution  parameters,  since  in  situ  calibration  will  probably  be 
necessary  for  accurate  waveform  studies. 
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Abstract 

The  Lopez  Island  OBS  Intercomparison  Experiment  provided  a  data 
set  of  sufficient  spatial  density  to  allow  study  of  the  propagation  of 
shot-generated  Stoneley  waves  as  well  as  ambient  background  noise. 

The  Stoneley  waves  were  observed  propagating  at  velocities  of  20  to 
50  m/sec.  Phase  velocities  were  determined  by  fitting  peaks  in  the 
frequency  wave  number  spectrum.  Croup  velocities  were  calculated 
by  narrowly  filtering  the  data  and  determining  the  arrival  time  of 
the  peak  in  the  frequency  packet.  Particle  displacement  plots 
illustrate  the  surface  wave  character  of  these  waves.  The  analysis 
of  the  ambient  background  noise  failed  to  produce  a  clearly  defined 
dispersion  curve  yet  it  did  allow  bounds  to  be  placed  on  the  velocities 
(20  to  50  m/sec).  The  data  were  modeled  using  eleven  layers  overlying 
a  half-space.  The  results  indicated  that  the  top  7  m  of  the  sediment 
column  at  Lopez  Island  is  best  approximated  by  two  zones.  In  the 
upper  zone  there  is  a  fairly  rapid  range  of  shear  velocity  with  depth. 
This  zone  overlays  a  region  in  which  the  shear  velocity  gradient  is 
much  lower.  Deep  ocean  background  noise  from  University  of  Washington 
ocean  bottom  seismometers  was  also  examined.  Although  insufficient 
data  precluded  any  velocity  analysis,  definite  similarities  exist 
between  these  data  and  noise  data  observed  at  Lopez  Island. 


i n t  rodnc  t ion 

The  Lopez  Island  OBS  Intercomparison  Experiment  provided  a  unique 
set  of  data  for  evaluating  the  particle  motion,  frequency  spectra,  and 
wave  number  spectra  of  background  seismic  noise  and  shot-generated, 
low-velocity  dispersive  surface  waves.  These  waves  have  been  named 
granddaddy  waves  (Heron  at  al.,1968)  because  of  their  extremely  large 
amplitude  relative  to  other  shot-generated  phases.  It  is  the  purpose 
of  this  paper  to  show  that  these  granddaddy  waves  and  the  Lopez  noise 
at  least  in  the  1-  to  5-Hz  band  propagate  as  Stoneley  waves  and  to 
suggest  that  deep  water  noise  in  this  frequency  band  may  also  propagate 
in  this  mode. 

Understanding  the  causes  and  modes  cf  propagation  of  ambient  deep 
ocean  bottom  noise  in  the  1  to  100  Hz  band  is  an  important  problem  to 
ocean  bottom  seismology  for  a  number  of  reasons.  First,  as  most  marine 
seismic  refraction  experiments  are  conducted  in  this  frequency  band, 
the  existence  of  any  significant  amount  of  noise  can  severely  degrade 
the  data.  Unpublished  OBS  data  suggest  that  noise  level  in  this  band 
varies  widely  with  geographic  location.  Second,  knowledge  of  noise 
characteristics  should  furnish  information  constraining  OBS  design 
parameters,  allowing  the  effects  of  noise  to  he  minimized. 

F’rior  to  the  Lopez  Island  Experiment  there  were  some  (unpublished) 
indications  that  deep  sea  noise  might  be  of  short  wavelength.  This 
being  the  case,  an  array  of  closelv  spaced  sensors  would  be  required 
to  determine  phase  velocity.  As  it  turned  out,  the  Lopez  Island  array, 
which  had  a  maximum  dimension  of  about  40  m,  was  reasonably  well-suited 
for  noise  analysis  (see  Sutton  et  al.,  1980,  for  locations).  As 
comparable  data  from  a  deep  water  array  do  not  now  exist,  only  the 
similarities  between  Lopez  Island  and  deep  water  data  car.  be  noted  and 
the  probable  characteristics  of  deep  water  noise  inferred. 

In  the  following  four  sections  the  nature  of  the  granddaddy  waves 
as  well  an  ambient  background  noise,  observed  both  at  Lopez  Island  and 
at  deep  ocean  sites,  will  be  addressed.  It  will  be  shown  that  the 
granddaddy  waves  observed  in  both  environments  have  similar  spectra 
and  velocity  dispersion  curves,  and  represent  Stoneley  waves 
propagating  along  the  sod iment-wu ter  interface.  Also,  it  will  be  shown 
that  noise  from  1  to  5  Hz  at  Lopez  Island  propagates  in  a  Stoneley 
mode,  which  will  lead  to  the  inference  of  similar  properties  for  deep 
sea  noise;  also  theoretical  models  will  be  compared  to  observed  data. 


Cha  rart eristics  of  Granddaddy  Waves  Produced  bv  Explosions 

Lopez  Island.  For  a  description  of  the  Lopez  Island  experiment  the 
reader  is  referred  to  Sutton  et  al.  0980).  During  this  experiment  a 
number  of  cap  shots  and  airgun  shots  were  fired  that  produced  low- 
veloeity  dispersed  wave  trains.  An  example  of  these  waves  is  shown  in 


_ 
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Figure  C-l  on  a  time-distance  graph.  Because  of  their  extremely  low 
velocity  relative  to  the  P  wave  velocity  of  water  and  the  proximity  of 
the  shots,  the  P  arrival  time  has  been  used  as  the  origin  time.  Ranges 
or  shot  locations  relative  to  the  array  were  determined  with  a  non-linear 
least  squares  algor  itlim,  which  solved  for  shot  location  and  velocity. 

The  inputs  for  this  technique  were  instrument  location  and  arrival  times 
for  a  given  frequency  packet.  For  a  more  complete  discussion  of  this 
technique,  see  Tuthill  (1980). 

The  dispersed  wave  trains,  as  illustrated  in  Figure  0-1,  have 
a  group  velocity  of  approximately  20  m/sec .  The  energy  in  these  waves 
is  concent  rated  in  a  1 -  to  5-!lz  band  as  shown  by  the  sample  spectrum 
in  Figure  0-2  from  airgun  shot  120.  The  calculated  Fourier  spectra 
from  instruments  with  1  -  Hr.  and  4.5-Hz  seismometers  show  little 
difference,  suggesting  that  the  low-frequency  cutoff  at  about  1  Hz 
may  be  real  and  not  entirely  caused  by  instrument  response.  It  was 
also  observed  that  these,  waves  were  highly  attenuated  on  a  hydrophone 
placed  about  1  m  above  the  bottom. 

Because  the  particle  motion  is  a  discriminant  o  f  the  wave  type,  a 
method  of  graphically  displaying  this  motion  was  devised.  Any  two 
components  are  combined  into  an  instantaneous  displacement  vector  and 
these  are  plotted  as  a  function  of  time.  Figures  (I- la  and  (1—  '3  b 
illustrate  the  synthesized  particle  motions  for  four  different  phase 
angles  and  three  different  amplitude  ratios,  between  two  components. 
Rectilinear  motion  corresponds  to  a  phase  difference  of  O'1  and 
elliptical  motion  corresponds  to  a  phase  difference  of  90°,  With  the 
actual  data,  ground  velocity  has  been  changed  t o  ground  displacement 
by  the  usual  frequency  domain  method.  An  example  of  these  aethods 
applied  to  the  dispersed  wave  train  generated  bv  airgun  shot  182 
is  shewn  in  Figures  (1-4  and  C-5.  This  plot  demonstrates  that  the 
particle  motion  is  elliptical  and  thus  is  consistent  with  the  motion 
of  a  surface  wave. 

In  determining  the  group  velocity  dispersion  curve  for  these  waves 
the  data  were  first  narrow-hand  filtered  (.5-Hz  hand  width)  at  various 
center  frequencies  with  a  4-polo  Butterworth  filter.  The  envelope  of 
this  filtered  wave  train  was  determined  bv  taking  the  modulus  of  the 
analytic  signal  defined  by  f ( t )  -  i  Fq(t),  where  FH(t)  is  the  Hilbert 
transform  of  f(t)  (Bracewell,  1965).  The  travel  time  for  this  frequency 
band  is  given  bv  arrival  of  the  peak  of  the  envelope.  An  example  of 
this  method  applied  to  the  spike  vertical  standard  for  the  capshot, 
decree  2'39,  is  shown  in  Figure  G-6.  These  times  furnished  the  travel 
time  input  for  the  non-linear  least  squares  method  described  above. 

The  determination  of  phase  velocity  proved  to  be  more  difficult. 

The  approach  was  to  take  the  frequency  wave  number  spectrum  of  the 
two-dimensional  array  and  use  the  peaks  in  the  wave  number  spectrum  for 
a  particular  frequency  to  establish  i.i/k  (the  phase  velocity).  In  this 
method  accurate  inter-element  spacings  are  critical  and  one  assumes 
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plane  waves  propagating  across  the  array.  In  actuality  the  shots  were 
close  enough  to  the  array  that  this  assumption  was  violated  signifi¬ 
cantly  and  the  method  did  not  yield  stable  results.  To  improve  the 
stability  the  array  element  positions  were  projected  onto  a  plane 
including  the  shot,  and  the  two-dimensional  frequency-wave  number 
spectrum  was  calculated.  The  resulting  w-k  values  were  contoured  and 
and  example  is  shown  in  Figure  G-7.  The  dashed  line  shows  the  peak 
that  was  fit  to  give  and  expression  for  o)(k).  The  final  observed 
group  and  phase  velocity  dispersion  curves  are  shown  in  Figure  G-8. 

Deep  Sea.  Observations  of  dispersive  Stoneley  waves  on  the 
sea  floor  have  been  made  several  times.  Hamilton  et  al.  (1969)  used 
submersibles  to  deploy  small  charges  at  short  distances  from  geophones 
to  measure  Stoneley  wave  velocities  and  hence  to  infer  sediment  shear 
velocities.  However,  the  experiment  described  by  Davies  (1965)  is 
more  comparable  to  the  results  discussed  here.  In  that  experiment  an 
array  of  hydrophones  was  deployed  on  the  sea  floor  to  record  shots 
detonated  on  the  sea  floor  at  ranges  up  to  a  few  kilometers.  Large- 
amplitude  dispersed  waves  were  recorded  that  Davis  interpreted  to  be 
Stonejey  waves.  In  this  case  the  use  of  hydrophones  to  record  these 
waves  was  successful  because  they  were  placed  on  the  sea  floor.  If 
they  has  been  positioned  a  few  meters  above  the  bottom  they  would 
probably  have  only  weakly  observed  these  waves.  An  example  of  time 
series  recorded  by  Davies  at  a  range  of  D.6  km  is  shown  in  Figure  G-9. 
The  spectrum  of  this  signal  is  displayed  in  figure  G-10  and  shows  a 
strong  peak  from  1  to  5  Hz.  (Note:  The  time  series  was  obtained  from 
photographically  enlarging  a  journal  figure  and  hence  the  spectrum 
obtained  from  this  time  series  is  only  approximate). 

The  resulting  group  velocity  dispersion  curve,  inferred  phase 
velocities,  and  model  that  fit  these  data  are  shown  in  Figure  G-ll. 
Note  that  the  dispersion  is  critically  dependent  on  the  shear  velocity 
gradient  in  the  upper  20  of  sediment. 


Characteristics  of  Ocean  Bottom  Noise 

Lopez  Island.  Numerous  samples  of  background  noise  were  taken 
during  the  Lopez  Island  experiment  (see  Sutton  et  al.,  1980).  Almost 
without  exception  the  noise  as  recorded  on  the  standard  instruments 
was  peaked  from  1  to  5  Hz.  An  example  from  the  plate  standard  is 
shown  in  Figure  G-12,  which  included  the  particle  velocity,  the  spectrum 
of  the  particle  velocity,  and  the  particle  displacement.  The  noise 
peaks  from  20  to  30  Hz  can  be  related  to  mechanical  sources  such  as 
motor  boats  and  are  not  included  in  this  analysis. 

To  establish  that  the  noise  is  propagating  in  the  Stoneley  wave 
mode  it  needs  to  be  verified  that  (1)  the  particle  motion  is  elliptical, 
(2)  the  amplitude  is  attenuated  in  the  water,  and  (3)  the  phase  velocity 
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has  the  correct  value  relative  to  the  sediment  shear  velocity.  The 
shear  velocity  of  the  sediment  can  be  established  from  the  shot¬ 
generated  dispersive  waves. 
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Particle  motion  diagrams  from  the  noise  displacement  as  recorded 
on  the  plate  are  shown  in  Figure  G-13.  This  figure  shows  clearly  that 
the  motion  is  predominantly  elliptical  for  the  2-Hz  frequencies  in  this 
record.  Pressure  recordings  of  the  noise  by  the  hydrophone  show  the 
1-  to  5-Hz  energy  to  be  highly  attenuated  relative  to  other  frequencies 
(see  Sutton  et  al.,  1980). 

The  phase  velocity  of  the  noise  as  a  function  of  frequency  was 
attempted  by  taking  the  three-dimensional  frequency  wave  number 
spectrum.  This  allows  one  to  determine  both  the  azimuthal  variation 
of  the  noise  and  its  velocity  at  a  particular  frequency.  Numerically 
the  wave  number  spectrum  is  very  sensitive  to  sensor  location  and  to 
the  assumption  of  plane  wave  propagation.  The  latter  assumption  would 
be  violated  if  the  noise  sources  are  close  to  or  within  the  array. 
Further,  the  wave  number  resolution  and  the  response  of  the  array  will 
not  be  ideal  (i.e.,  like  a  delta  function.  Fig.  G-lAa).  Some  or  all 
of  these  factors  are  responsible  for  the  degradation  of  the  actual 
wave  number  spectra  obtained.  An  example  of  the  spectrum  taken  at 
2  Hz  is  shown  in  Figure  G-14b.  However,  in  spite  of  the  poor 
resolution,  we  can  definitely  bound  the  phase  velocity  at  2  Hz  between 
20  and  50  m/sec. 


Deep  Sea  Noise 

Unfortunately,  data  adequate  to  define  the  velocities  of  deep  sea 
noise,  that  is,  data  having  closely  spaced  sensors,  do  not  now  exist. 

The  best  that  can  be  done  at  present  is  to  document  the  spectra  and 
particle  motion.  To  do  this  we  have  chosen  a  worst  case  sample  of 
noise  taken  at  Lat.  16°N,  Long.  145°W.  This  was  the  site  of  a  joint 
experiment  in  which  both  Scripps  and  the  University  of  Washington 
deployed  several  instruments.  All  the  OBS's  at  this  general  location 
were  placed  on  sediment  and  had  very  high  levels  of  noise  at  1  to  3  Hz. 
In  fact  the  noise  was  so  large  that  in  order  to  extract  refraction 
arrivals  the  data  had  to  be  severely  filtered  to  remove  this  noise 
component.  This  filtering,  of  course,  severely  degraded  the  information 
of  interest.  An  example  of  the  ground  noise  velocity,  its  spectrum, 
and  corresponding  ground  displacement  are  shown  in  Figure  G-15.  These 
data  show  a  very  high  noise  peak  at  about  2Hz.  The  corresponding 
particle  motion  in  Figure  G-16  shows  the  noise  to  be  predominantly 
elliptical  in  motion  and  non-uniform  in  azimuth.  Although  not  proof 
that  this  noise  is  propagating  as  Stoneley  waves,  this  is  suggested 
as  a  possibility  based  on  the  similarities  with  the  Lopez  noise  and 
the  shot- generated  dispersed  waves. 
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For  purposes  of  comparison,  the  deep  sea  noise  spectrum  is  compared 
with  the  other  spectra  in  Figure  C-17.  The  similarities  are  quite 
obvious. 


Comparison  with  Theory 

In  this  section  it  is  shown  that  the  data  presented  in  this  paper 
are  consistent  with  the  propagation  of  a  dispersed  surface  wave  at 
the  interface  between  the  water  and  the  sea  floor.  Further,  it  is 
shown  that  the  dispersion  characteristics  are  related  to  the  gradient 
in  sediment  shear  velocity  just  below  the  sea  floor. 

A  liquid  layer  over  an  elastic, homogeneous  half-space.  A  convenient 
starting  point  is  with  the  simple  model  in  order  to  evaluate  the  effect 
of  the  liquid  layer  on  the  velocity  and  dispersion.  Biot  (1952) 
investigated  this  problem  and  showed  that  Stoneley  waves  were  simply 
the  high-frequency  limit  of  Rayleigh  waves,  or  equivalently,  for 
wavelengths  short  compared  to  the  water  thickness  the  fundamental  mode 
surface  waves  are  confined  to  the  water-solid  interface.  The  equations 
relating  the  phase  velocity  Vp  to  the  water  thickness  H  and  the  wave¬ 
length  ''  are 


4(1  -  (Vp/p?>‘V-(2  -  (Vp/fiP)V/(l  "  (V/a2) = 

ll _ IVlT. _  tan[2"h((V  /nj)2  -  1)J] 

o2  (<Vp/«i>-  -  l)1  X  P 


and 


for  V  >  a 
P 


4(1  -  (V n/R;)2)5-  ((2  -  V  /B?)?)r/((1  -  V  /n?)?): 


P  (V  /Bt)1*  „  ’ 

_L  P  ?  tan  I.l2*h  ((V  /a , ) -  1)") 

P  2  ((Vp/rM)?  -  I)' 

for  V  <  i, 

P  1 


where  r*.^ ,  ot2  are  compressional  wave  velocities  of  the  water  and  solid, 
&2  Is  the  shear  velocity  of  the  solid,  pj  and  P2  are  the  densities,  and 
X  is  the  wavelength. 


Deep  Sea 
Stoneley  Waves 


ectra  from  Lopez  shot,  Davies  shot,  Lopez  noise,  and  deep  sea  noise. 


These  equations  have  been  solved  for  the  fundamental  mode  for 
various  values  of  83’  h  and  1,  and  the  results  are  listed  in  Table  G-l . 
It  is  clear  from  the  table  that  as  32  approaches  ctj  the  waves  are 
highly  dispersed  whereas  for  82  the  waves  are  essentially  non- 

dispersive.  It  is  also  apparent  that  the  phase  velocities  are 
insensitive  to  water  depth  for  82  Physically  this  is  due  to 

the  wavelengths  in  the  solid  being  so  small  relative  to  the  water 
thickness  that  the  wave  is  effectively  confined  to  the  interface. 

It  should  he  noted  that  for  the  frequencies  of  interest  (1  to  10  Hz) 
and  for  82  =  3U  m/sec  even  a  1-m  water  layer  is  unimportant,  and  these 
waves  are  non-di spersive.  The  non-dispersive  nature  is  due  to  the 
homogeneous  half-space  assumption  and  implies  that  the  observed 
dispersion  must  be  caused  by  velocity  changes  below  the  sediment  - 
water  interface.  This  is  the  same  conclusion  reached  by  Davies  (1965). 
The  effect  of  the  water  in  the  high-frequency  limit  is  to  depress  the 
phase  velocities  to  about  0.86  83  compared  to  0.95  83  for  Rayleigh 
waves  on  the  free  surface  of  an  incompressible  half-space. 

A  layered  half-space.  For  the  case  of  a  single  water  layer  over- 
lying  a  homogeneous  half-space  we  have  shown  that  for  82  «1  the 

water  layer  can  be  effectively  disregarded,  as  far  as  dispersion  is 
concerned.  This  implies  that  the  observed  dispersion  at  Lopez  Island 
and  in  the  deep  ocean  must  be  caused  by  velocity  changes  in  the  top 
f<w  meters  of  the  sediment.  This  is  how  Davies  (1965)  modeled  his  data 
and  we  have  undertaken  .1  similar  approach  for  the  Lopez  Island  data. 

A  program  to  calculate  Rayleigh  wave  group  velocity  for  a  layered 
half-space  was  used  to  generate  a  model  that  fit  the  observed  group 
velocities  reasonably  well.  The  final  layered  model  (Fig.  G-18)  can 
be  divided  into  two  zones  of  nearly  equal  thickness  ('^3  m).  The 
upper  zone  has  a  high  shear  velocity  gradient  and  the  lower  zone 
has  a  very  low  shear  velocity  gradient. 

We  note  that  this  general  model  is  the  same  type  used  by  Davies 
(1965)  although  the  sediment  shear  velocities  inferred  from  Lopez 
Island  are  appreciably  lower.  At  Lopez  Island  the  shear  velocity  in 
the  upper  few  meters  are  about  15  to  20  m/sec.  There  was,  however,  no 
Airy  phase  seen  in  the  data,  so  further  modeling  is  warranted. 


Gone  bis  i on_s 

We  have  shown  that  large-amplitude,  low-velocity  dispersed  waves 
generated  bv  shots  at  Lopez  Island  correspond  to  Stoneley  wave 
propagation  at  the  water-sediment  interface  and  that  the  dispersion  is 
caused  by  the  shear  velocitv  structure  in  the  upper  few  meters  of  the 
sediment.  We  find  that  the  sediment  shear  velocities  at  Lopez  Island 
are  about  a  factor  of  two  less  than  those  found  by  Davies  in  a  deep 
water  sample.  Clearly  the  shear  velocity  will  be  determined  by  the 


c;- 29 


sediment  type,  porosity,  and  degree  of  lithification,  which  in  general 
will  depend  on  the  geographic  location. 

We  have  also  shown  that  the  background  noise  at  Lopez  Island  has 
spectra,  particle  motion,  and  phase  velocities  similar  to  the  shot- 
generated  dispersed  waves,  which  implies  that  the  predominant  noise 
is  also  propagating  as  Stoneley  waves.  In  deep  water  the  documentation 
ot  noise  is  insufficiently  precise  to  prove  that  Stonely  waves  are  the 
predominant  mode  of  propagation,  but  the  available  data  are  not 
inconsistent  with  this  hypothesis. 


There  are  several  implications  of  these  findings  that  are  of 
importance  to  ocean  bottom  seismology. 

(1)  The  very  low  velocity  of  the  noise  implies  wavelengths  of 
a  few  meters  to  tens  ot  meters  depending  on  the  actual  shear  velocity 
of  the  sediments.  Consequent ly, the  noise  will  be  incoherent  between 
sensors  spaced  more  than  a  few  tens  of  meters  apart. 


(2)  The  fact  that  the  noise  is  most  likely  propagating  in  the 
Stonele'7  mode  implies  that  the  amplitude  of  the  noise  will  decrease 
as  e^”“'  away  from  the  interface  (Z  =  distance  from  the 

interface,  '  =  wavelength').  For  frequencies  of  2  Hz  and  velocities  of 
20  m/sec  e'  “  f  v  =  e  '  for  7.  "V.  1.6  m.  This  implies  that  the 
noise  can  be  greatly  reduced  either  by  placing  a  hydrophone  a  few 
meters  above  the  bottom  or  burying  a  seisnu  meter  a  few  meters  below 
the  bottom. 


Another  aspect  of  these  data  that  we  have  not  yet  discussed  is 
the  narrow  spectrum  of  the  dispersed  waves,  the  energy  falling 
between  7  and  about  5  Hz.  We  suspect  that  the  high-frequency  cutoff 
may  he  controlled  by  non-elastic  attenuation  in  the  sediments  and  the 
low-frequency  cutoff  controlled  by  the  shear  velocity  profile  in 
the  sediments.  For  low  frequencies,  or  long  wavelengths,  the  waves 
may  be  sampling  a  higher  shear  velocity  which  will  cause  them  to  be 
decoupled  from  the  main  dispersed  wave  train.  These  combined  effects 
may  concentrate  the  noise  power  from  2  to  5  Hz.  Since  most  OBS's  are 
placed  in  sediment  ponds  and  most  long-range  refraction  profiles  use 
this  frequency  hand  for  propagation  of  the  "signal",  it  is  clear 
that  improvement  in  the  design  of  OBS's  could  result  in  a  significant 
increase  of  signal-to-noise  ratio. 
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Abs  tract 


Data  from  the  1978  Lopez  Island  OBS  Intercomparison  Experiment  and 
deep  sea  data  from  University  of  Washington  OBS's  show  that  there  is  a 
considerable  amount  of  waveform  distortion  resulting  from  the  conversion 
of  horizontal  motion  into  vertical  motion,  here  called  cross-coupling 
distortion.  This  distortion,  which  substantially  reduces  the  significance 
of  waveform  matching  with  synthetic  seismograms,  appears  to  result  from 
rotation  imparted  to  the  OBS  package  by  near-vertically  traveling  shear 
energy.  The  degree  of  this  rotation  seems  to  depend  on  the  instrument 
surface  area  above  the  seafloor  and  the  geometry  and  surface  area  of  the 
"feet"  connecting  the  package  to  the  seafloor.  The  sensitivity  and 
response  of  the  seismometers  within  the  package  to  this  rotation  depends 
on  the  precise  location  of  the  seismometers  with  respect  to  the  axis  of 
rotation.  The  results  suggest  how  to  modify  OBS  designs  to  minimize  these 
effects . 


H-2 


Introduction 


One  of  the  principal  objectives  of  the  Lopez  Island  OBS  Intercomparison 
Experiment  was  an  evaluation  of  the  coupling  of  the  various  OBS  designs  to 
the  ocean  bottom.  As  a  result  of  this  experiment  it  was  found  that  most 
instruments  (except  the  standards)  had  a  cross-coupling  problem  in  which 
near-vertically  traveling  shear  motion  was  coupled  into  vertical  motion. 

This  effect  was  also  observed  during  the  transient  tests,  in  which 
horizontal  transients  of  force  produced  substantial  outputs  on  the 
vertical  seismometers,  and  it  is  documented  by  Sutton  et  al .  (1980).  In 
this  paper  we  further  document  this  effect  by  using  controlled  sources  at 
the  Lopez  site  and  in  deep  water  sites,  and  in  addition  we  offer  possible 
explanations  and  remedies  for  its  cause. 


Cross-Coup!  ing  _Da_t_a _ f  rom  Lope z  Island  Shots 

One  of  the  best  demonstrations  of  the  cross-coupling  problem  is  the 
data  from  a.irgun  shot  182  (Fig.  H-l).  There  are  three  distinct  energy 
packets  on  these  shot  records.  Before  1  sec  on  Figure  H-l  there  is  a 
P  arrival  distinguished  by  its  predominant  vertical  motion.  Between  1  sec 
and  5  sec  the  motion  is  predominantly  horizontal,  due  to  near-vertically 
traveling  S  waves.  The  beginning  of  this  wave  train  is  probably  leaky 
S  waves.  Beyond  a  travel  time  of  6  sec  the  motion  is  due  to  a  dispersive 
surface  wave  whose  characteristics  are  described  by  Tuthill  et  al .  (1980). 
It  is  the  S  part  of  the  record  in  which  we  are  interested  here,  as  recorded 
on  the  vertical  components.  OBS  recordings  of  this  shot  are  shown  in 
Figure  H-2.  We  note  that  whereas  the  standard  instrument's  component 
showed  very  little  motion  for  the  S  phase,  the  OBS's  show  a  wide  range  of 
vertical  response  to  this  wave  train.  This  response  is  clearly  an  artifact 
of  the  instrumentation  and  is  a  distortion  of  the  true  particle  motion. 

A  relative  ranking  of  the  instruments  according  to  the  degree  of  distortion 
is  shown  in  Figure  H-3  by  plotting  the  ratio  of  the  maximum  vertical 
amplitude  of  the  leaky  mode  to  the  amplitude  of  the  first  two  cycles  of 
the  P  wave. 


Cross-Coupling  Data  from  Deep  Ocean  Results 

In  tli is  section  we  document  a  case  of  deep  water  cross-coupling 
distortion  similar  to  that  seen  in  the  Lopez  Island  data.  The  deep  water 
data  come  from  a  deployment  of  four  University  of  Washington  OBS's  on  the 
Cascadia  Basin  off  the  coast  of  Washington  in  September  1978.  The  water 
depth  at  that  site  was  about  3  km,  the  sediment  thickness  is  about  2  km, 
and  the  sediments  are  regular  and  flat-lying.  The  OBS's  were  deployed  in 
a  linear  array  about  4  km  long;  three  of  the  OBS's  had  concrete  anchors 
and  one  had  a  tripod  anchor.  Schematic  diagrams  of  these  instruments  are 
shown  in  Figures  H-4  and  H-5.  Vertical  component  records  for  one  shot 
across  the  array  are  shown  in  Figure  H-6  and  horizontal  component  records 
for  the  same  shot  are  shown  in  Figure  H-7.  Figure  H-8  shows  a  comparison 
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Fig.  H-2.  Recordings  of  airgun  shot  182  by  the  OBS's  and  the 
standard  verticals.  The  large  amplitude  energy  between  2  and 
3  sec  seen  on  most  of  the  OBS’s  is  horizontal  motion  cross 
coupled  in  the  vertical  seismometers. 


Fig.  H-3.  Ranking  o£  Che  instruments  according  to  the  degree  of  cross  coupling.  The 
ordinate  represents  the  ratio  of  the  maximum  amplitude  of  the  cross-coupled  shear 
motion  to  the  first  two  cycles  of  P  motion. 
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Fig.  H-5.  A  schematic  drawing  of  (he  University  of  Washington  OBS  with  the 
"tripod"  anchor.  CB  is  the  center  of  buoyancy;  CM  is  the  center  of  mass 
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Fig.  H-7.  A  comparison  of  horizontal  component  recordings  of  the  same  shot  shown  in 
Figure  H-6.  The  largest  amplitude  events  correspond  to  refracted  shear  arrivals. 
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of  vertical  and  horizontal  data  for  a  concrete  anchor  and  the  tripod 
anchor  for  this  shot.  From  an  inspection  of  Figures  H-6,  H-7,  and  H-8,  it 
is  clear  that  the  P  arrival  has  a  high  coherence  across  the  array  and  is 
essentially  independent  of  the  anchor  type.  However,  when  we  inspect  the 
refracted  S  part  of  the  wave  train  we  note  that  although  the  horizontals 
are  reasonably  coherent  the  vertical  motion  shows  little  correlation  from 
OBS  to  OBS.  This  is  especially  seen  when  comparing  the  concrete  and 
tripod  anchors.  We  interpret  this  to  be  waveform  distortion  due  to 
cross  coupling  of  horizontal  motion  into  vertical  motion,  in  a  manner 
exactly  analogous  to  the  Lopez  Island  results. 

We  point  out  that  we  have  chosen  this  data  set  to  demonstrate  this 
result  because  the  S  and  P-S  arrivals  are  well  separated  from  the  P 
arrival  because  of  the  large  sediment  thickness.  In  cases  of  much  thinner 
sediments  (a  few  meters  to  tens  of  meters),  the  P-S  arrivals  will  occur 
during  the  P  arrivals  and  the  cross  coupling  will  considerably  distort 
the  P  waveform,  making  it  uninterpretable. 


Possible  Causes  of  Cross-Coupling  Distortion 

In  Figure  H-3  we  ranked  the  (IBS's  according  to  degree  of  cross  coupling 
and  we  note  that  those  instruments  with  a  large  cross  section  in  the  water 
and  a  low  spread  of  support  on  the  ocean  bottom  (inverted  pendulums)  have 
the  worst  degree  of  cross-coupling.  Tripod-tvpe  instruments  with  large 
foot  pads  spread  far  apart  do  a  lot  better,  whereas  the  spike,  plate  and 
instruments  with  remote  sensor  packages  do  the  best.  This  in  itself 
suggests  that  the  cause  of  the  cross  coupling  is  rocking  induced  by  lateral 
translation  caused  by  shear  motion. 

There  are  several  mechanisms  by  which  rocking  motion  can  be  translated 
into  vertical  motion.  Three  mechanisms  are  described  below  and  in  each  of 
these  we  assume  that  the  hydrodynamic  resistance  caused  by  the  cross  section 
in  the  water  is  greater  than  the  spring  constant  associated  with  the 
coupling  to  the  soft  sediments.  That  is,  we  assume  that  the  cross  section 
in  the  water  remains  fixed  with  respect  to  the  water  and  the  lateral  trans¬ 
lation  is  taken  up  by  depressing  the  feet  into  the  sediment.  The  three 
mechanisms  are: 

(a)  Fnifrom  sediment  properties  under  each  foot.  In  this  case  the 
rocking  will  occur  about  an  axis  close  to  the  center  of  area  above  water 
and  the  effect  on  the  vertical  and  horizontal  seismometers  will  depend  on 
the  precise  location  of  the  sensors  with  respect  to  the  axis.  This  is 
denoted  in  Figure!!-^  as  motion  Rj  and  has  been  experimentally  verified  by 
rocking  the  Fniversitv  of  Washington  tripod  OBS  about  axes  at  different 
heights  from  the  vertical  seismometer.  Depending  on  the  location  of  the 
vertical  seismometer,  frequencies  at  twice  the  driving  frequency  can  be 
induced;  this  mode  is  especially  important  when  the  sensor  is  at  the  axis 
of  rotation. 
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POSTULATED  CROSS  COUPLING 
MECHANISMS 


SHEAR  MOTION 


Fig.  H-9.  A  schematic  diagram  of  possible  cross-coupling  mechanisms. 
Rj  represents  rotation  caused  by  rocking  on  a  uniform  bottom. 

R.,  represents  rocking  caused  by  a  non-uniform  bottom,  represented 
by  different  spring  constants  S]  and  S ? .  The  sensitivity  of  the 
vertical  seismometer  to  R,  depends  on  its  precise  location  with 
respect  to  the  axis  of  rotation. 
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(b)  Non-uniform  sediment  properties  under  each  foot.  This  case  is 
equivalent  to  having  different  spring  constants  under  each  foot  and  can  be 
caused  by  the  instrument  landing  unevenly  on  a  bottom  with  a  high  shear 
strengh  gradient  in  the  upper  few  centimeters.  In  this  case,  rocking  will 
be  effectively  about  the  foot  with  the  highest  spring  constant  and  is 
denoted  by  motion  R2  in  Figure  H-9.  Again  the  effect  will  depend  on  the 
precise  location  of  the  sensor  but  in  general  it  will  produce  dominantly 
vertical  motion. 

(c)  Inverted  pendulum  case.  These  instruments  usually  have  the 
sensors  close  to  the  mud  line  and  the  center  of  buoyancy  and  area  a  meter 
or  so  above  the  mud  line.  Mechanisms  (a)  and  (b)  are  not  applicable  to 
this  case  since  it  is  impossible  with  small  rotations  about  a  center  of 
rotation  near  the  flotation  to  produce  vertical  to  horizontal  cross¬ 
coupling  ratios  of  one  or  greater,  as  has  been  observed.  The  most  likely 
cross-coupling  mechanism  in  this  case  is  a  dynamic  interaction  of 
buoyancy  forces  with  the  spring  constant  associated  with  bottom  coupling. 
More  specifically,  small  rotations  caused  by  lateral  translation  of  the 
foot  will  induce  buoyant  restoring  force  acting  against  the  bottom 
coupling  spring  constant.  Resonances  at  the  exciting  frequency  could 
produce  large  amplification  of  the  vertical  motion.  We  have  not  attempted 
a  mathematical  or  physical  modeling  of  this  situation.  To  be  complete  such 
a  model  should  incorporate  hydrodynamic  effects. 


Con  c  1  u  s  i  o  n_s 

We  have  demonstrated  that  instrumental  waveform  distortion  of 
near-vert ica 1 ly  traveling  shear  energy  occurs  on  all  DBS's  to  some  degree. 
The  distortion  results  from  the  cross-coupling  of  horizontal  into  vertical 
motion  and  is  dependent  on  the  geometry  of  the  instrument,  the  location  of 
the  sensors  within  the  instrument,  and  the  homogeneity  of  the  ocean  bottom. 

The  cross  coupling  can  he  reduced  by  having  small  sensor  packages 
remote  from  the  main  instrument  or  by  having  low-profile  instruments  with 
widely  spaced  feet  and  large  footpads.  These  factors  should  be  considered 
along  with  the  other  DBS  design  parameters  described  by  Sutton  et  al. 

(1980)  for  coupling  frequencies  and  Tuthill  et  al.(1980)  for  noise  reduction. 
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Abstract 

First  arrival  compressional  wave  signals  from  an  airgun  source,  as 
detected  bv  a  varietv  of  seismometers  in  a  shallow  bay,  are  remarkably 
uniform;  however,  minor  variations  in  wavelet  appearance  imply  some  com¬ 
bination  of  instrument  response  and  coupling  to  the  bottom.  Signal 
spectra  show  typically  a  spectral  peak  at  12  Hz  and  an  envelope  very 
similar  to  that  expected  from  an  airgun  source.  Those  instruments  with 
a  decoupled  geophone  package  have  spectra  most  like  the  theoretical 
spectrum  but  spectra  for  the  other  instruments  are  not  significantly 
different.  Little  variation  exists  in  spectra  between  tripod-mounted  and 
inverted-pendulum  OBS  configurations  for  the  low-amplitude  P-waves 
observed  here.  The  signal  source  is  the  principal  influence  on  the  re¬ 
sulting  spectra  rather  than  OBS  configuration  or  bottom  coupling. 
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Introduction 


The  Lopez  Island  intercomparison  experiment  provides  a  unique  set 
of  data  to  study  bottom  coupling  and  response  of  ocean  bottom  seismographs 
(OBS)  to  seismic  waves  in  shallow  water.  The  data  are  particularly  useful 
because  signals  from  each  instrument  were  brought  ashore  in  real  time  by 
cable  and  recorded  on  a  24-channel  digital  recorder  to  provide  wide- 
dynamic-range  records  of  multiple  instruments  for  direct  comparison  during 
extended  computer  analysis  in  the  laboratory.  One  of  the  seismic  sources 
recorded  during  the  experiment  was  a  40-inch^  bolt  airgun  fired  to  the 
OBS  array  from  various  azimuths  and  ranges.  Because  an  OBS  is  commonly 
used  as  a  receiver  for  an  airgun  source  during  marine  experiments  it  is 
useful  to  examine,  under  controlled  conditions,  the  effectiveness  of  such 
an  instrument  for  the  detection  of  airgun-generated  ground  motion  at  a 
water-solid  interface.  The  interaction  of  seismic  waves  in  shallow  water 
and  an  instrument's  response  on  the  bottom  is  complex  (Lewis  and  Tuthill, 
1980;  Tuthill  et  al.,  1980)  and  in  this  report  we  will  restrict  our 
attention  to  compressional  waves  generated  at  ranges  of  300  to  700  m. 


Data 


Figure  1-1  shows  a  set  of  airgun  shots  from  a  seismic  profile  in  a 
NW  direction  across  Shoal  Bay  on  Lopez  Island  as  detected  by  the  well- 
coupled  spike  vertical  standard  geophone  and  by  the  Oregon  State  University 
OBS  vertical  component.  Two  wave  groups  occur  within  the  first  two  seconds 
of  record  and  have  phase  velocities  of  550  and  2,850  m/s.  The  shots  in 
Figure  1-1  are  about  80  m  apart  at  ranges  between  300  to  700  m  from  the 
OBS  array.  The  two  record  sections  are  nearly  identical  except  for  phase 
and  amplitude  differences  (instrument  responses  have  not  been  equalized). 

In  order  to  quantize  characteristics  of  the  seismic  signals  from  the 
various  instruments,  we  have  analyzed  the  first  arrivals  from  two  shots 
at  700-  and  600-m  distance.  It  was  necessary  to  use  two  shots  because 
not  all  instruments  were  wired  to  the  digitizer  at  the  same  time.  We 
assert  that  the  proximity  of  the  shot  points  to  each  other  and  the  uni¬ 
formity  of  the  bottom  (Sutton  et  al. ,  1980b)  result  in  comparable  seismic 
signals.  We  restrict  our  attention  to  the  high-velocity  arrivals  since 
the  lower-velocity  arrivals  may  represent  a  shear  phase.  At  the  greatest 
distance  commensurate  with  adequate  signal  level,  we  study  the  first 
1.024  s  (256  samples)  of  signal  from  each  of  the  instruments  as  shown  in 
Figure  1-2.  Spike,  neutral  density,  and  plate-mounted  vertical  geophone 
standards  (SFV,  NDV,  PLV)  represent  the  range  of  strong  to  weak  bottom 
coupling.  The  three  or  four  letter  code  in  the  signal  label  identif  es 
the  originating  institution  for  each  instrument. 


Theoretical Airgun  Signal 

The  analysis  that  follows  requires  a  source  signal  from  the  airgun 
used  during  the  experiment.  Because  such  a  signal  was  not  recorded  during 
tiie  experiment,  we  generated  a  theoretical  signal  by  using  the  method  of 


Record  sections  of  arrivals  to  the  spike  vertical  standard  and  to  the  vertical  component 
State  OBS  along  line  6,  which  was  shot  northwest  from  the  array.  Shot  125  is  approximate. 
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-W*^V\A/'WW'W\/^  — w  NDV  AG  125 


PLV  AG  125 
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Ziolkowski  (1970).  The  theory  is  based  on  dynamic  and  thermodynamic 
relations  derived  by  Gilmore  (1952)  for  an  oscillating  spherical  bubble. 

The  theoretical  damped  pressure  wave  generated  under  conditions  at  Lopez 
Island  is  given  in  Figure  l-3a.  Reflections  of  the  signal  from  the  air- 
water  and  water-bottom  interfaces  were  delayed,  attenuated  and  summed 
until  their  amplitude  became  less  than  1%  of  the  initial  value  (assuming 
1/r  spreading)  as  shown  in  Figure  l-3b.  The  signal  was  filtered  at  125  Hz 
and  resampled  at  4  ms  (250  sps)  to  duplicate  the  conditions  under  which  the 
OBS  signals  were  digitized  (Fig.  I-3c).  A  fast-Fourier  transform  (FFT) 
algorithm  generated  the  power  spectrum  of  Figure  3d  from  1.024  s  (256 
samples)  of  the  filtered  airgun  signal.  The  spectrum  has  a  peak  at  14.5  Hz 
corresponding  to  the  69-ms  bubble  pulse  period  shown  in  Figure  I-3a. 


Observed  Spectra 

The  initial  portion  of  the  observed  wave  signals  in  Figure  1-2  were 
processed  by  the  FFT  to  generate  the  spectra  shown  in  Figures  1-4,  1-5 
and  1-6.  Figure  1-4  shows  the  theoretical  spectrum  to  be  in  good  agree¬ 
ment  with  the  spectra  for  the  standard  geophones.  In  general,  the 
observed  spectra  have  a  single  maximum  at  about  13  Hz  and,  neglecting 
minima  due  to  interference  effects,  fall  off  very  much  like  the  theoreti¬ 
cal  spectra.  For  comparison.  Figure  1—4  shows  a  shot  close  to  the  spike 
vertical,  which  is  nearly  identical  to  the  theoretical.  A  geophone  on 
land  recorded  shot  188  also,  but  the  high  background  noise  tends  to  mask 
any  spectral  peak. 

Spectra  for  eight  instruments  are  shown  in  Figure  1-5  for  airgun 
shot  125.  The  narrowest  peaks  are  displayed  by  MITV  and  HIGBV ,  both  of 
which  had  geophones  in  a  separate  pressure  case  to  isolate  the  sensor 
from  recording  vibration.  The  low  profile  of  the  fall-away  package 
appears  to  have  the  added  advantage  of  producing  more  faithful  representa¬ 
tion  of  the  ground  motion  in  comparison  to  the  spectra  of  other  OBS  s. 

Tli is  is  probably  due  to  the  improved  coupling  of  the  sensor  case  of  the 
bottom  (Sutton  et  al.,  1980a).  The  responses  for  the  remaining  five 
instruments  are  shown  in  Figure  1-6  for  airgun  shot  192.  These  spectra 
are  noisier  than  those  for  shot  125  but  demonstrate  a  peak  corresponding 
to  the  theoretical  signal.  The  USGS  instrument  shows  a  lack  of  energy 
above  10  Hz.  The  reason  for  it  is  unknown,  but  the  waveform  is  markedly 
unlike  the  others. 


Discussion 

The  waveforms  as  detected  by  the  OBS  instruments  for  shots  125  and 
192  are  in  general  similar  to,  but  in  detail  somewhat  different  from, 
those  of  the  standard  instruments.  These  differences  are  probably  due  to 
instrument  responses  that  have  not  been  incorporated  in  the  presentation 
of  seismograms.  Comparison  of  spectra  from  distant  shots  (which  should 
minimize  cross  coupling  between  horizontal  and  vertical  motions  of  the 
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Fig.  1-6.  Seismograms  and  spectra  for  five  instruments 
at  a  distance  of  600  m.  The  spike  vertical  standard 
for  comparison  with  SPV  AG125  in  Figure  4. 


from  airgun  shot  192 
is  shown  in  upper  left 
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instrument  duo  to  rocking)  are  similar  and  resemble  the  spectrum  expected 
on  a  theoretical  basis.  However,  subtle  differences  remain  between 
instrument  spectra.  'Hie  spectra  most  consistent  with  the  theoretical  are 
those  whose  geophone  packages  coupled  properly  to  the  bottom,  i.e.,  the 
external  geophone  package.  Of  the  remaining  OBS  types,  there  is  little  in 
the  spectra  of  the  compress ional  waves  to  distinguish  signals  of  tripod- 
mounted  spheres  from  those  of  inverted  pendulum  designs  (where  the  float 
is  located  higher  than  the  sensor).  Exceptions  to  this  are  the  UCSB 
instrument  whose  instrument  response  peaks  at  8  Hz  and  the  USGA  instru¬ 
ment  mentioned  earlier. 

Particularly  noticeable  in  '-'iguro  1-5  is  the  frequency  of  the 
spectral  peak  at  12  Hz  for  the  observed  signals  in  comparison  to  14.5  Hz 
for  the  theoretical  signal.  If  we  assume  that  this  is  due  to  frequency- 
dependent  attenuation  of  the  signal,  we  can  calculate  the  Q  of  the  medium. 
The  spectral  ratio  between  the  theoretical  signal  and  the  spike  vertical 
was  formed  for  airgun  shot  125.  The  values  of  this  ratio  versus  frequency 
are  shown  plotted  in  Figure  1-7  and  a  best  fit  straight  line  is  drawn 
through  the  points.  The  slope  of  the  line  is  related  to  attenuation  and 
we  compute  a  Q  of  8.1  for  the  2.85  km/s  seismic  layer. 

The  basement  rock  of  Lopez  Island  is  a  siltstone-sandstone 
conglomerate  of  marine  origin  (Sutton  et  al . ,  1980b).  The  value  of  Q 
agrees  favorably  with  laboratory  measurements  of  Toksoz  et  al.  (1979)  on 
Boise  and  Navajo  sandstone  (•>  -  6.9  -  7.3).  We  conclude  that  the  frequency 
shift  of  the  spectral  peaks  at  Lopez  Island  is  due  to  the  low  Q  of  the 
basement  rock. 

These  observations  of  air gun  signals  confirm  studies  of  OBS  response 
from  lift  tests  (Sutton  et  al.,  1980a)  and  from  particle  motion  studies 
(Lewis  and  Tuthill,  1980)  but  are  somewhat  loss  definitive.  For 
low-amplitude  P  waves  having  a  predominately  vertical  ray  path,  seismograms 
obtained  with  a  wide  variety  of  OBS  designs  compare  favorably  with  signals 
from  the  standards.  Large-amplitude  waves  with  considerable  horizontal 
motion  adverselv  affect  the  OBS  signal  except  for  OBS  with  externally 
deployed  geophones.  P>roadening  and  splitting  of  airgun  spectra  may  be 
caused  by  coupling  and  the  physical  configuration  of  the  instrument.  The 
spectra  observed  for  low-amplitude  P-waves  suggest  that  the  spectrum  of 
the  source  dominates  over  either  OBS  configuration  or  coupling  of  he 
instrument  to  the  bottom.  Finally,  we  suggest  that  the  spectra  of  first 
arrival  P-waves  can  be  determined  adequately  from  the  signals  of  most 
of  the  OBS  instruments  at  Lopez  Island. 
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Abstract 

High-amplitude,  narrow  band  noise  that  correlates  with  periods  of 
high  ocean  bottom  currents  and  the  tidal  cycle  is  occasionally  observed 
on  ocean  bottom  seismometers  (OBS).  The  geophones  on  OBS's  of  different 
configurations  are  not  equally  sensitive  to  this  noise  and  hydrophones 
are  almost  unaffected.  With  a  suitable  design,  it  should  be  possible  to 
eliminate  this  noise  problem. 


Introduction 


The  effects  of  water  currents  on  DBS's  have  been  discussed  in  the 
literature,  most  recently  by  Kasahara  et  al.  (in  press).  They  relate  the 
observed  vibrations  to  {Carman  vortex  shedding  off  various  parts  of  the 
instrument  and  find  that  effects  of  Karman  shedding  are  observable  at 
water  speeds  above  about  10  cm/sec.  This  value,  however,  is  probably 
dependent  on  configuration  of  the  package.  In  Karman  shedding  vortices  of 
water  spin  away  downstream  from  an  obstruction  and  exert  a  force  on  the 
object.  The  vortices  spin  off  at  well-defined  intervals  resulting  in  a 
harmonic  force  on  the  body  .  The  frequency  of  the  force  is  governed  by 
the  equation  f=cj,  where  f  is  the  frequency,  is  the  speed,  d  is  the 
diameter  of  the  body  and  C  is  the  strouhal  number.  The  strouhal  number 
depends  on  the  shape  of  the  body  and  varies  between  0.15  and  0.2  in  this 
case.  The  shedding  is  non-linear  in  that  the  vortices  do  not  begin  to 
spin  off  until  a  critical  speed  is  reached.  The  motions  resulting  from 
these  forces  can  be  amplified  if  the  shedding  frequency  is  near  a  natural 
resonance  of  the  obstruction.  Karman  shedding  in  air  caused  by  wind 
blowing  around  wires  or  flag  poles  generates  a  moan  or  whistle.  Water 
density  is  about  1000  times  greater  than  that  of  air;  thus  dynamic 
pressures  against  a  body  in  water  are  about  1000  times  greater  than  for 
that  body  in  air  at  the  same  speed.  The  pressure  exerted  by  a  3-m/sec  (6  - 
kt)  wind  is  approximately  equivalent  to  a  current  of  10  cm/sec.  since 
pressure  increases  as  the  square  of  the  speed. 

In  this  paper  we  describe  the  effects  of  current  on  HIG  OBS's, 
discuss  design  parameters  that  should  minimize  this  noise,  and  speculate 
about  the  effects  of  currents  on  long-period  seismic  observations. 


Analysis 

The  correlation  of  noise  level  recorded  on  an  OBS  south  of  the  island 
of  Hawaii  in  5  km  of  water  with  current  speed  and  tide  can  be  seen  in 
Figure  J-l.  The  noise  level  has  been  reduced  to  nm  (rms)  of  equivalent 
vertical  ground  motion  to  show  the  sensitivity  of  this  OBS  to  currents. 

The  tide  is  the  theoretical  ocean  tide  (Longman,  1959)  computed  for  a 
point  about  halfway  between  the  OBS  and  the  current  meter  (which  were 
separated  by  about  50  km).  The  lack  of  a  perfect  correlation  between 
noise  level  and  current  speed  is  probably  the  result  of  two  factors:  (1) 
the  distance  separating  the  two  instruments  and  (2)  tae  sensitivity  of  the 
OBS  to  currents  from  different  directions.  The  OBS  is  elongate  in  one 
direction  (Fig.  J-2)  and  may  be  more  sensitive  to  current-induced  noise  in 
one  direction  than  in  the  other.  Spectra  of  the  noise  observed  by  the 
horizontal  geophone  in  this  OBS  are  shown  in  Figure  J-3 .  The  lower 
spectrum  was  taken  when  no  current  noise  was  visible  and  the  upper  when  it 
was  strong.  The  monochromatic  character  (with  overtones)  and  narrow 
bandwidth  of  the  current  noise  is  obvious.  The  frequency  of  the  noise  is 
observed  to  increase  slightly  as  the  noise  level  increases.  These  spectra 
are  similar  to  those  observed  by  Bradner  et  al.  (1965),  who  suggested 
that  the  noise  was  caused  by  Karman  vortices  from  a  radio  antenna. 
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We  first  connected  the  2-Hz  “whistle"  heard  on  HIG  DBS's  with  the 
tidal  cycle  in  1976,  and  also  noticed  that  the  resonant  frequency  of  our 
recovery  radio  beacon  antennas  was  close  to  2  Hz.  In  1977,  seven  OBS's 
were  dropped  with  various  antenna  configurations;  some  with  balls  attached 
near  their  centers  or  tops,  some  with  flags  attached  to  the  top,  and  some 
with  2-inch-diameter  PVC  pipes  around  them.  The  recordings  obtained  from 
these  OBS's  showed  some  changes  in  resonant  frequency  but  no  consistent 
differences  in  current  noise  susceptibility.  Although  the  antennas  may  be 
the  main  source  of  this  noise,  other  parts  of  the  OBS  package  may  also  be 
sources.  In  1978  we  designed  a  new  style  OBS  (Fig.  J-4)  with  isolated 
geophones  in  an  attempt  to  overcome  the  current  noise  problem,  to  isolate 
the  geophones  from  tape  recorder  noise,  and  to  improve  mechanical  coupling 
to  the  ocean  bottom. 

The  new  configuration  is  much  less  sensitive  to  current  noise. 

Figure  J-5  compares  noise  levels  observed  by  the  two  styles  of  HIG  OBS's 
located  within  100  m  of  each  other  in  deep  water.  The  spikes  on  each 
record  are  signals  from  explosives  set  off  for  a  refraction  line.  The 
strong  increase  in  the  background  level  on  the  right  side  of  traces  B  and 
C  is  current  noise  identical  to  that  described  above.  Note  that  the 
hydrophone  (trace  A)  and  the  two  lower  traces  (D  and  E)  are  not  affected 
by  this  noise.  The  hydrophone  is  attached  to  the  package  containing  the 
geophones  that  produced  traces  B  and  C.  The  configuration  of  the  OBS  that 
recorded  traces  D  and  E  (Fig.  J-4)  differs  significantly  from  the  OBS  that 
produced  traces  B  and  C  (Fig.  J-2)  in  that  the  geophones  are  mechanically 
decoupled  from  the  main  package.  The  geophom  package  thus  has  a  much 
lower  profile,  less  complexity,  a:  <i  smaller  cross  section  exposed  to  the 
current.  We  do  not  know  if  all  (or  any)  of  l he  above  differences  account 
for  the  reduction  in  current  no;  st:  sensitivity,  although  the  analysis  of 
Kasahara  et  al.  (in  press)  suggests  that  antennas  attached  to  the  sensor 
package  are  a  major  problem.  He  stales  that  the  problem  can  be  reduced  by 
attaching  a  fin  to  the  antenna. 

We  estimate  that  from  5*  to  i  52  ol  our  geophono  data  from  OBS's  of 
the  type  shown  in  Figure  J-2  is  degraded  by  current  noise.  While  the 
isola ted-geophone  package  also  has  other  advantages  (Sutton  et  al., 

1980a),  the  improvement  in  data  finality  shown  in  Figure  J-5  is  worth  the 
added  mechanical  complexity  of  the  isolated-geophone  design. 

Although  many  different  OBS  configurations  are  in  use,  only  a  few 
authors  have  mentioned  current  noise  problems  in  their  data.  Whether  it 
does  not  occur,  is  not  recognized  o:  observed,  or  is  simply  ignored  is 
important  to  OBS  design.  It  OBS's,  especially  OBS's  with  internal 
geophones,  now  exist  that  are  insensitive  to  current  noise,  then  one  of 
the  motives  for  building  isolated-gcophone  OBS's  becomes  unimportant.  One 
of  the  purposes  of  the  Lopez  Island  OBS  Intercomparison  Experiment  was  to 
determine  the  effects  of  current  noise  on  the  various  instruments  (Sutton 
et  al.,  1980b).  Unfortunately,  the  maximum  currents  observed  during  the 
experiment  (5  cm/sec)  were  less  than  the  10  cm/sec  found  by  Kasahara  et 
al.  (in  press)  as  a  threshold  for  the  onset  of  K arm an  vortex  shedding. 

The  noise  level  at  I.opez  during  the  quietest  times  was  about  50  nm  at  2 
Hz,  which  is  higher  than  most  of  the  current  noise  amplitudes  from  Figure 
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Fig.  J-4.  HIG  external  geophone  OBS .  The  geophones  are  deployed 
away  from  the  main  recording  package  to  decrease  current  and 
tape  recorder  noise  and  to  improve  coupling  to  the  bottom. 
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Fiji.  I-j.  One-minute  average  noise  levels  verses  time  for  the  two 
different  styles  of  OBS  reeorded  in  the  same  place  at  the  same 
time  during  the  shooting  of  a  refraction  line.  Trace  A  is  the 
record  frorr.  a  hydrophone;'  B  is  f rom  a  vertical;  and  C  is  from 
a  horizontal  geophone  in  an  one  1 osed-geophone  OBS  (Figure  J-2)  . 
n  and  K  are  records  f rom  vertical  and  horizontal  geophones  in 
a  deployed-go. -phone  OBS  (Figure  J-4) .  The  increase  in  noise 
level  to  the  right  of  traces  B  and  C  is  caused  by  bottom  currents. 
This  noise  does  not  appreciably  affect  the  hydrophone  (A)  or  the 
deployed  geophones  (D  and  F.)  . 
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J-l*  Even  if  noise  was  generated,  it  would  be  barely  detectable.  An  HIG 
OBS  identical  to  that  shown  in  Figure  J-2,  and  known  to  be  susceptible  to 
current  noise,  was  present  at  the  Lopez  Experiment  and  recorded  none.  The 
Lopez  Experiment  did  not  test  OBS's  for  current  noise  susceptability  nor 
has  any  test  for  this  parameter  been  made  for  most  of  the  OBS's  now  in 
use . 


Vortex  shedding  is  not  the  only  possible  mode  of  coupling  current 
noises  to  OBS's.  An  instrument  resting  on  soft  sediments  can  be  tilted  by 
the  current  pressure  against  it.  For  example,  an  OBS  with  a  circular 
cross  section  about  one  meter  in  diameter  centered  about  0.5  m  above  the 
ocean  bottom  with  a  tripod  base  and  mass  of  200  kg  in  Lopez-type  sediments 
(Sutton  et  al.,  1980b)  would  tip  about  5  x  10-6  m  (10“5  radius)  in 
response  to  a  10  cm/sec  current.  While  a  tilt  such  as  this  would  not  be 
important  for  short-period  seismometers,  long-period  horizontal 
seismometers  would  be  driven  off  scale.  A  small  fluctuation  in  current 
speed  would  produce  objectionable  noise  at  shorter  periods. 


Conclusions 


Current  noise  is  a  known  problem  for  at  least  some  ocean  bottom 
seismometers.  At  short  periods,  the  principal  mode  of  noise  transmission 
appears  to  be  the  shedding  of  Karman  vortices  from  antennas  and  other 
resonant  obstructions.  At  long  periods,  tilting  of  the  OBS  caused  by 
ocean  bottom  currents  could  be  a  serious  problem.  Both  problems  can  be 
improved  significantly  by  isolating  the  sensors  from  possible  current 
noise  sources  (analogous  to  planting  geophones  on  land  away  from  trees) 
and  by  keeping  the  cross  section  to  the  water  as  small  as  possible.  The 
isolated-sensor  HIG  OBS  appears  to  have  a  relatively  low  current  noise 
sensitivity.  Most  OBS's,  however,  have  never  been  tested  for  current 
noise  sensitivity. 
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ground  motion  much  as  a  low  pass  filter  does,  and  high  buoyancy  tends  to 
counteract  this  effect.  These  effects  are  observed  in  the  Lopez  data,  which 
consist  of  signals  from  mechanical  transient  tests,  cap  shots,  airgun  pulses, 
and  general  background  noise.  The  modification  of  response  is  pronounced  for 
some  instruments  and  barely  noticeable  in  others.  Instruments  that  stand 
high  in  the  water  relative  to  their  base  width  tend  to  be  susceptible  to 
rocking  motion  that  shows  up  as  a  mechanical  cross  coupling  between  horizontal 
and  vertical  motion.  Correlation  of  Lopez  results  with  coupling  theory 
suggests  that  i t  is  possible  to  design  ocean  bottom  seismometers  that  will 
couple  well  to  any  sediment.  Current  levels  at  the  Lopez  site  (<5  cm/ sec) 
wore  too  small  to  produce  noticeable  effects  on  any  of  the  instruments; 
however,  the  same  design  criteria  that  will  minimize  coupling  problems  will 
also  lessen  problems  caused  by  ocean  currents. 
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